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; ABSTRACT 
J"""-. , We present spectral energy distributions (SEDs), Spitzer colours, and infrared (IR) luminosi- 

ty"} ■ ties for 850 /im selected galaxies in the Great Observatories Origins Deep Survey (GOODS) 

' Northern field. Using the deep Spitzer Legacy images and new data and reductions of the 

, VLA-HDF radio data, we find statistically secure counterparts for 60 per cent (21/35) of our 

submm sample, and identify tentative counterparts for another 12 objects. This is the largest 
sample of submm galaxies with statistically secure counterparts detected in the radio and with 
Spitzer. Half of the secure counterparts have spectroscopic redshifts while the other half have 
photometric redshifts. We find that in most cases the 850 /im emission is dominated by a sin- 
JL \ gle 24 /im source, with a median flux density of 241 /iJy, leading to a median 24 /im to 850 /im 

flux density ratio of 0.040. A composite rest-frame SED shows that the submm sources peak 
' at longer wavelengths than those of local ultraluminous infrared galaxies (ULIRGs). Using 

, a basic greybody model, 850 /im selected galaxies appear to be cooler than local ULIRGs 

L* ■ of the same luminosity. This demonstrates the strong selection effects, both locally and at 

. | high redshift, which may lead to an incomplete census of the ULIRG population. The SEDs 

, of submm galaxies are also different from those of their high redshift neighbours, the near-IR 

5_i ■ selected BzK galaxies, whose mid-IR to radio SEDs are more like those of local ULIRGs. Us- 

ed 1 ing 24 /im. 850 /im and 1.4 GHz observations, we fit templates that span the mid-IR through 

radio to derived the integrated IR luminosity of the submm galaxies and find a median value 
of Lir(8 — 1000/im) = 6.0 x 10 12 Lq. By themselves, 24 /im and radio fluxes are able to pre- 
dict Lin reasonably well because they are relatively insensitive to temperature. However, the 
submm flux by itself consistently overpredicts Ltr when using spectral templates which obey 
the local ULIRG temperature-luminosity relation. The shorter Spitzer wavelengths sample the 
stellar bump at the redshifts of the submm sources, and we find that the Spitzer photometry 
alone provides a model independent estimate of the redshift, er( Az/(1 + z)) = 0.07. The me- 
dian redshift for our secure submm counterparts is 2.0. Using X-ray and mid-IR data, only 5 
per cent of our secure counterparts (1/21) show strong evidence for an active galactic nucleus 
(AGN) dominating the IR luminosity. 

Key words: galaxies: formation - galaxies: evolution - galaxies: starburst - infrared: galaxies 
- submillimetre 



1 INTRODUCTION 

Deep extragalactic surveys at X-ray through radio wavelengths 
have revealed different galaxy populations at high redshift. Sev- 
eral of these galaxy populations, including Lyman Break Galax- 



ies (LBGs, Steidel et al. 1998), massive BzK galaxies (Daddi et 
al. 2004) and bright submillimetre (submm) galaxies (see Chap- 
man et al. 2005 and references therein), are thought to contribute 
significantly to the universal star formation at z ~ 2 — 3. Whether 
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these and other high redshift populations are completely different 
galaxy systems or simply galaxies seen at different phases in their 
evolution requires a complete multi-wavelength picture. 

Since the LBG and BzK galaxy populations are selected in 
high resolution optical images, determining their counterparts at 
other wavelengths is straightforward, and so they have been studied 
extensively in other wavebands (e.g. Shapley et al. 2005; Daddi et 
al. 2005). However, galaxies detected with the Submillimetre Com- 
mon User Bolometer Array (SCUBA, Holland et al. 1999) on the 
James Clerk Maxwell Telescope (JCMT) suffer from a large beam 
size (~15 arcseconds at 850 fim). The density of optical galaxies in 
deep Hubble Space Telescope (HST) surveys is such that there are 
often ~10 optical galaxies in every SCUBA beam, making identi- 
fication difficult. 

Follow-up studies of submm galaxies present two main chal- 
lenges: 1) obtaining deep enough data at other wavelengths to have 
confidence that the plausible counterpart(s) has been detected; and 
2) deciding among possible counterparts in order to make the cor- 
rect identification. The first of these requires considerable amounts 
of telescope time, while the second can be particularly difficult 
given the large beam size. Deep radio observations have proven to 
be the best way to localize the submm emission coming from high 
redshift 850 pm selected galaxies (e.g. Ivison et al. 2000; Smail 
et al. 2000; Barger, Cowie, & Richards 2000). Sources detected in 
deep 1.4 GHz VLA observations are rare enough that the proba- 
bility that they are randomly associated with a submm galaxy is 
quite low (Ivison et al. 2002; Borys et al. 2004b). However, even in 
the deepest radio maps, there is still a substantial fraction (at least 
1/4) of submm galaxies which remain undetected, and 10 per cent 
have more than one radio source to choose from (Ivison et al. 2002; 
Chapman et al. 2005). Our current knowledge of submm galaxies is 
biased toward the radio-detected sub-sample which appear to con- 
sist of massive objects (Borys et al. 2005; Greve et al. 2005) at 
z ~ 2 — 3 (Chapman et al. 2005). These galaxies are rare and have 
very high star formation rates (Lilly et al. 1999; Scott et al. 2002), 
and hence may represent an early phase in the evolution of massive 
elliptical galaxies. 

To complement, or as an alternative, to the radio, Multiband 
Imaging Photometer for Spitzer (MIPS; Rieke et al. 2004) obser- 
vations with the Spitzer Space Telescope can probe the rest-frame 
mid-IR emission in submm galaxies. In addition to being much 
closer in wavelength to the submm than the optical, the surface 
density of 24 /im detected sources is such that the probability of 
randomly finding a source within a SCUBA search radius is much 
less than in the optical. 

The road to determining the contribution of submm galaxies 
to the global star formation at high redshift requires a detour to ac- 
curately estimate the IR luminosity for individual galaxies. In order 
to get a measure of the total IR luminosity and star formation rate 
for high redshift galaxies, an extrapolation of the flux at near- or 
mid-IR wavelengths is often required. Another avenue to get to the 
total infrared luminosity is through the radio flux and the radio- 
infrared correlation, which is well established locally but has not 
been fully tested at high redshift (e.g. Condon 1992 and references 
therein; Appleton et al. 2004). Both of these approaches are indi- 
rect and require assumptions about the shape of the SED which is 
known to differ for different galaxy types (e.g. Dale et al. 2005). 
Before Spitzer, the region of the IR SED shortward of 850 pm was 
unexplored in high redshift submm galaxies (Blain et al. 2002). In 
order to obtain the most accurate measure of the IR luminosity re- 
quires at least several data points from 8-1000 /im restframe, and 
we can now begin to do this with Spitzer. It will not be until the 



Herschel Space Observatory (Pilbratt 2001) that we are able to fill 
in the most crutial gap in the far-IR SED with deep, high resolution 
photometry at ~ 200/im. 

In this paper, we use the deep Spitzer images from the GOODS 
Legacy programme and a new reduction of the 1.4 GHz VLA ra- 
dio data (Morrison et al. in preparation) to study a large sample of 
bright (Ssbo > 2 mjy) submm galaxies in the GOODS-N field. We 
present SEDs, Spitzer colours, and IR luminosities for these galax- 
ies. We explore how different points on the IR SED probe the IR 
luminosity and use the stellar SED to investigate the presence of 
AGN in these systems. 

The format of this paper is as follows. Section 2 describes 
the submm sample in GOODS-N. Section 3 summarizes the multi- 
wavelength data in the field and describes the new Spitzer observa- 
tions. A summary of our multi-wavelength idenfications is given in 
Section 4, including statistics to assess the reliability of the associ- 
ations. In Section 5 we briefly describe the SED templates used in 
this paper. Section 6 presents the Spitzer properties of this submm 
sample and discusses how these properties evolve with redshift, 
while Section 7 presents a composite IR to radio SED for our sam- 
ple. Infrared luminosities, as estimated from the 24 /im, 850 pm 
and 1.4 GHz flux densities, as well as star formation rates, are pre- 
sented in Section 8. Section 9 looks briefly at whether there is a 
significant contribution from an AGN to the IR luminosity of our 
submm sources. We summarize our main points in Section 10. The 
appendix contains tables of positions and flux densities for all the 
submm sources and their counterparts, notes on individual objects 
and postage stamp figures at optical, mid-IR and radio wavelengths. 

All magnitudes in this paper use the AB system unless oth- 
erwise noted. We assume a standard cosmology with Ho = 
72kms _1 Mpc -1 , Q M = 0.3 and Ov = 0.7. 



2 SUBMILLIMETRE SAMPLE 

The SCUBA 'super-map' of GOODS-N contains 40 submm 
sources detected above 3. 5a at 850 pm (Borys et al. 2003, Pope et 
al. 2005). Table Al lists the complete sample of GOODS-N submm 
mapping sources, along with the submm positions, the 850 pm 
flux densities and the coordinates of the identified galaxy counter- 
parts (see Sections 3-6). Since all mapping observations taken in 
GOODS-N are combined into the super-map, the depths achieved 
in this map are very non-uniform. The raw noise levels of the de- 
tected sources range from 0.3-4.1 mjy, as seen in Table Al. The 
submm observations and data analysis for this sample are described 
in detail in Borys et al. (2003, 2004b) and Pope et al. (2005). For 
simplicity we will refer to these three papers as Paper I, Paper II and 
Paper III, respectively. We use the naming convention with the pre- 
fix 'GN' as listed in Paper III when referring to individual sources 
in this paper, however, the full submm names are listed in Table 
Al. The SCUBA data that we use here are identical to those used 
in Paper III, within which the image and noise map can be found. 

Since we have chosen a 3.5a detection threshold, we have to 
be aware of the possibility of spurious sources in our sample. Due 
to the effects of confusion and noise, sources detected in SCUBA 
maps at low significance have typically had their fluxes boosted by 
some factor which depends on both the source flux and local noise 
level (see Coppin et al. 2005 and references therein). Several ways 
of dealing with this bias have been used in other submm studies; the 
fundamental issues one needs to deal with are that low signal-to- 
noise sources are, on average, at somewhat lower flux than they ap- 
pear in maps, and that sources with higher noise (at a given signal- 
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to-noise level) are more likely to be spurious. In order to obtain the 
most robust sample, we choose to apply the flux deboosting pre- 
scription of Coppin et al. (2005) which is a Bayesian method de- 
pendent on prior knowledge of the source count model. Applying 
this prescription 1 , we find that five of our 40 sources have a non- 
negligible probability (> 5 per cent) of deboosting to zero flux. 
The threshold chosen is based on having a very low false detection 
rate in our final catalogue. As expected (Ivison et al. 2002; Cop- 
pin et al. 2005), these five sources, GN27, GN29, GN33, GN36 
and GN38, have the highest noise levels. Given that we expect 
some small fraction of sources in our sample to be spurious, we 
excluded these five from our list in order to make a more robust 
sample. Therefore the final submm sample in GOODS-N contains 
35 850 fim sources. We note that the main results of this paper are 
not strongly affected by details of how we construct our final sam- 
ple. Three of the five sources which are excluded from the sam- 
ple based on flux deboosting (GN27, GN29 and GN38) appear to 
have radio and/or Spitzer identifications and are probably genuine 
submm sources. We list these additional sources at the bottom of 
Table Al. However, we stress that they are not in our main cata- 
logue of 35 submm sources, the selection of which is based solely 
on the SCUBA data, uninfluenced by whether there appears to be a 
counterpart at other wavelengths. 

The flux deboosting algorithm also provides corrected values 
of the flux and noise for all remaining sources, based on the proba- 
bility distribution of flux. Obtaining the most accurate submm pho- 
tometry is important in order to correctly model the SED and inves- 
tigate variations in properties with submm flux. For the 35 sources 
in our sample we list the raw and deboosted fluxes in Table I All It 
is clear that the low signal-to-noise ratio and high noise sources are 
most affected by flux-boosting. 

None of these 35 sources is individually detected above 3cr in 
the SCUBA 450 fim map. However, the sample is statistically de- 
tected at 6.3 ± 2.2 mJy. This value is obtained by measuring the 
450 fim flux at the positions of the 850 fim detections and calcu- 
lating a noise-weighted average. We expect that this measurement 
will be biased low due to positional uncertainties (from low signal- 
to-noise ratios) in both the 850 /im and 450 fim maps. This value is 
low compared to what is predicted from SED templates (roughly by 
a factor of 3). In addition to the positional offsets, the 450 fim cal- 
ibration is very uncertain and the noise properties in deep 450 fim 
SCUBA maps are known to be complicated. Hence stacking flux 
down to these levels is essentially untested at 450 fim. For this rea- 
son (and the reasons discussed in Paper III, where we abandoned 
the attempt to detect individual sources) we do not consider 450 fim 
data further. 



3 GOODS MULTI- WAVELENGTH DATA 

The GOODS survey consists of three major components: Chan- 
dra 1 Msec X-ray observations (Alexander et al. 2003b); deep HST 
optical imaging using the Advanced Camera for Surveys (ACS) 
in four bands (Giavalisco et al. 2004); and deep Spitzer imaging 
at five infrared wavelengths (Dickinson et al., in preparation). In 
addition to the extensive space-based imaging campaign, several 
ground-based programmes are also targeting the GOODS fields for 

1 We used a simplified approach of adopting a prior distribution coming 
from a single chop function. This does not strictly apply to the entire super- 
map (where several chop strategies have been used), but we verified that 
this makes little difference to the results (see Coppin et al. in preparation). 



imaging and spectroscopy. Optical counterparts for a large fraction 
of our submm sample using the HST images of GOODS-N (Gi- 
avalisco et al. 2004) are presented in Paper III. The Spitzer imaging 
of GOODS-N was completed in November 2004 and since it is the 
focus of much of the analysis in this paper, we discuss the observa- 
tions and data analysis in more detail here. 

3.1 Spitzer observations and data 

Spitzer observations of the field were obtained as part of the 
GOODS Legacy programme (Dickinson et al., in preparation). The 
entire GOODS-N field was imaged at 3.6, 4.5, 5.8 and 8.0 /im with 
InfraRed Array Camera (IRAC; Fazio et al. 2004) and at 24 pm 
with MIPS. These observations are currently the deepest Spitzer 
images. 

The Spitzer/GOODS data reduction is presented in Dickinson 
et al. (in preparation). The IRAC data have a resolution of ~ 2 arc- 
seconds. It was found that using a 'Mexican hat' kernel for IRAC 
source detection improved deblending in crowded regions. This oc- 
casionally resulted in improved IRAC source positions and fluxes 
for the submm counterparts. The IRAC photometry was measured 
using matched apertures in SExtractor (Bertin & Arnouts 1996). 
For the present analysis, we use a 4 arcsecond diameter aperture 
and apply aperture corrections as determined through simulations 
of the GOODS IRAC data. The uncertainties in the IRAC pho- 
tometry were also estimated through these simulations, since the 
SExtractor values were not accurate. The formal la point source 
sensitivities of the IRAC data for isolated sources are 0.026, 0.044, 
0.290 and 0.321 ^Jy at 3.6, 4.5, 5.8 and 8.0 ^m, respectively. Based 
on simulations, the 50 per cent detection completeness limits for the 
SExtractor IRAC catalog used in this paper are 0.4 ji]y at 3.6 and 
4.5 fim and 0.9 fily at 5.8 and 8.0 ^m. The IRAC and ACS (vl.0) 
images are aligned to an overall RMS of 0.25 arcsec and we have 
corrected for the known offset of -0.38 arcsec in declination be- 
tween the ACS and radio frames 2 . All positions in this paper are 
relative to the radio frame. 

The MIPS 24 fim data were reduced using an IDL pipeline, 
and photometry was obtained by fitting the Point Spread Func- 
tion (PSF). MIPS source detection was carried out using the IRAC 
source positions as a prior, which facilitates the comparison be- 
tween the IRAC and MIPS catalogues. 24 fim sources are selected 
to be above 3a. Since the 24 catalogue is fundamentally de- 
fined by the existence of IRAC priors, the low signal-to-noise ratio 
24 fim sources are much more likely to be real. Final astrometric 
accuracy is better than 0.2 arcseconds. More details on the MIPS 
source catalogues are given in Chary et al. (in preparation). The fi- 
nal 24 fim image reaches a la depth of about 5 fiiy, with an 84 per 
cent completeness limit of 24 fiSy. The number densities of sources 
in the Spitzer catalogues are listed in TableQ 

3.2 Spectroscopic and photometric redshifts 

Photometric redshifts have been estimated for a large sample of 
sources in GOODS-North using all available optical/near-IR data 
(U KPNO, B, V,R,I,z SUBARU, Capak et al. 2004, J, K s KPNO, 
-B435, V606, *775, 2850 ACS HST, Giavalisco et al. 2004). These 
extensive photometric data have been fit to a suite of SED templates 
to estimated redshifts (see Mobasher et al. 2004 for more details). 
In Paper III, we looked in detail at the photometry and redshift 

2 http://data.spitzer.caltech.edu/popular/goods/Documents/goods_drl.html 
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Table 1. Number density of sources in GOODS-N multi-wavelength cata- 
logues. The submm number densities are from the number counts given in 
Borys et al. (2003). 



Catalogue 


nfarcmin 2 ) 


Submm S > 6 mjy 


0.14 


Submm S > 2 mjy 


0.53 


Radio 4a 


1.8 


Radio 3a 


11.3 


MIPS 3a 


12.5 


IRAC 3cr 


88.9 


Radio and MIPS 3a 


2.2 



probability distribution for each submm counterpart to check for 
noisy or inconsistent values and only reported photometric redshifts 
which were well constrained. 

In addition to photometric redshifts for a large fraction of opti- 
cal galaxies, there are roughly 1500 spectroscopic redshifts over the 
160 square arcminute GOODS-N field (Cohen et al. 2000; Cowie et 
al. 2004; Wirth et al. 2004; Chapman et al. 2005). This means that 
several counterpart candidates already have spectroscopic redshifts 
(SectionlOl. 



3.3 New radio data reduction 

The GOODS-N region was initially imaged with the VLA at 8.5 
and 1.4GHz (Richards et al. 1998; Richards 2000). The original 
1.4GHz catalogue contained 371 objects (n ~ 0.8arcmin -2 ) and 
the RMS in the centre of the image was around 7.5 /iJy. The radio 
data used in this paper are a combination of the reprocessed VLA 
A-array HDF data plus new VLA B-array data (~28 hr, Morrison 
et al. in preparation). The reprocessed A-array data lead to a 25 per 
cent reduction in the phase centre, and the A+B array data lead to a 
30 per cent reduction in the phase centre noise, a ~ 5.3 /iJy. In ad- 
dition, the local RMS is significantly improved at greater distances 
from the centre of the map. Catalogues were made down to 4a and 
3a, and the number density of sources in each is given in TableQ 
The radio positional error is 0.16 and 0.26 arcsec for a 5a and 3a 
detection, respectively. These new data and reductions of the VLA- 
HDF radio data will be presented in Morrison et al. in preparation. 

In order to estimate the radio fluxes, we used three different 
resolution radio maps: the full resolution (1.6 arcsec) map; a 3 
arcsec convolved map; and a 6 arcsec convolved map. The lower 
resolution maps are useful for searching for lower surface bright- 
ness sources and for accurately measuring fluxes of the resolved 
sources (Owen et al. 2005). Of these three measurement, the full 
resolution map yielded the best signal-to-noise ratio in most cases. 
For sources that had a nearby radio source, we only used the full 
resolution map. We utilized the Astronomical Imaging Processing 
System (AIPS) JMFIT task to measure the flux density of the ra- 
dio sources, which are corrected for primary beam attenuation and 
bandwidth smearing effects. Sources which were found by JMFIT 
to have zero as the minimum size of the major axis were assumed 
to be unresolved. For such unresolved sources the fitted peak flux 
density for the Gaussian functions fitted with JMF I T is the best es- 
timate for the total flux density (Owen et al. 2005). Errors in the flux 
density and position were calculated using the method discussed in 
Condon (1997). 



4 MULTI-WAVELENGTH IDENTIFICATIONS 

To identify counterparts to our sample, we have used all available 
multi-wavelength data, in particular the new Spitzer images and the 
new reduction of the VLA 1 .4 GHz radio data. We employ a search 
radius of 8 arcsecond, which we derived by minimizing the prob- 
ability that K or more submm galaxies (out of M ) have at least 
one radio source within the search region (see Paper II for more de- 
tails). Note that we have slightly expanded the search radius from 
the 7 arcseconds used in Paper II and Paper III in light of the new 
radio reduction; when carrying out this calculation with the new 
radio catalogue, there is a clear minimum at 8 arcseconds. 

The counterpart identification of submm galaxies comes down 
to a set of statistical estimates. For all possible counterparts within 
the search region we need to determine the reliability of the asso- 
ciation. In the submm, this is usually done using Poisson statistics 
to calculate the probability of finding the object at random at that 
position (e.g. Lilly et al. 1999; Ivison et al. 2002). The probability 
that a source of a given flux density, S, is randomly found within 6 
of a submm source is given by 

PS = 1 - exp(-7rn (>s) 6> 2 ), (1) 

where n>s is the surface density of sources above flux density level 

5 per unit solid angle, and 9 is the search radius. The final proba- 
bility of random association usually includes a correction factor to 
take into account the fact that one is looking for counterparts within 
a given search region using catalogues which have a specific depth 
(Downes et al. 1986; Dunlop et al. 1989). 

In this paper, we denote the corrected Poisson probability (us- 
ing the method of Downes et al. 1986) as P. Since we expect a large 
fraction of the submm galaxies to be detected in the radio and mid- 
IR, we want to use as much information as possible when assessing 
the reliability of the associations. We made a joint catalogue by 
matching radio and MIPS sources within 1 arcsecond. Recall that 
since IRAC priors were used in making the MIPS catalogue, we al- 
ready have a joint MIPS and IRAC catalogue. For all counterparts 
discussed in this paper, we have calculated P using several cata- 
logues and criteria: 1) a > 4a radio catalogue (R); 2) a > 3a MIPS 
24 /im catalogue (M); 3) a joint radio and 24 /im catalogue (R/M); 
and 4) a joint 24 /im and red IRAC catalogue (M/I). In the latter, 
we use the MIPS/IRAC catalogue and make cuts on colour, rather 
than magnitude or flux, when calculating n in Equation^ We have 
assigned counterparts in a systematic way, starting with the cata- 
logue with the smallest number density and progressing to other 
catalogues if potential counterparts are not found. TableQlists the 
number densities for each of these catalogues at the limit of the sur- 
vey. Note that we confident in using these low signal-to-noise ratio 
catalogues since all of our potential counterparts are also detected 
at much higher significance in all four IRAC channels, therefore it 
is unlikely that they are random noise peaks in the MIPS or radio 
images. 

A counterpart is considered secure if it has P < 0.05. This is 
the same threshold used by Ivison et al. (2002; 2005) and ensures 
that the fraction of incorrect counterparts remains low. We also list 
tentative counterparts for sources which have 0.05 < P <0.20. Sta- 
tistically we expect only ~ 10 per cent of these to be just random 
associations, but clearly we cannot tell which ones are genuine. 
Therefore we have taken the conservative approach to only use the 
secure counterparts in the analysis in this paper, but note there are 
no large changes in the results if we included the tentative counter- 
parts. 

As an additional test into the reality of the secure counterparts, 
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we have simulated random positions in GOODS-N and applied the 
above procedure to look for counterparts with P < 0.05. The re- 
sults of these Monte Carlo simulations is that we find a counterpart 
with P < 0.05 about five per cent of the time therefore we expect 
that only one of our secure counterparts is a random association. 

Using the above procedure, we have found secure counterparts 
for 21 (out of 35) submm galaxies and an additional 12 tentative 
counterparts. Table lA2l provides the optical-through-radio photom- 
etry for all of the 33 identified submm galaxies, and the probability 
that the counterpart is a random association, P. All but one of the 
secure counterparts are detected in the radio, at 24 /im, and in all 
four IRAC channels, but this is partially an artifact of our identifi- 
cation procedure, since we require P < 0.05 in one of these multi- 
wavelength catalogues in order for a idenfication to be classified as 
secure. However, the presence of bright 24 /im or red IRAC coun- 
terparts enhances our confidence in low signal-to-noise ratio radio 
detections as submm counterparts. There are two sources for which 
we are unable to assign a unique identification due to the fact that 
they have multiple counterparts which are equally likely. 

The list of securely identified sources in GOODS-N includes 
GN14, also known as HDF850.1 (Hughes et al. 1998), which is 
discussed in detail in Dunlop et al. (2004, see also Appendix Al 
for more details). The submm position is known from deep MER- 
LIN/VLA data. However, because the submm source lies behind a 
bright elliptical galaxy, we cannot separate the Spitzer flux density 
of the submm counterpart from that of the elliptical. Lensing from 
the foreground elliptical galaxy may also complicate this submm 
system. Since we can only give upper limits to the flux density at 
mid-IR wavelenghts, we exclude it from the rest of the analysis 
in this paper. In addition, the Poisson probabilities favour another 
counterpart for GN14 (see Appendix Al), therefore GN14 demon- 
strates that a strict statistical approach can fail occasionally, and it 
is important to consider all available multi-wavelength data when 
assigning counterparts to submm galaxies. 

Ours is not the first study to combine SCUBA and Spitzer 
data and we now summarize what is currently known. Egami et 
al. (2004) looked at the coincidence of IRAC and MIPS detections 
with 10 submm sources from the 8-mJy SCUBA Survey (Scott 
et al. 2002; Fox et al. 2002) using the Lockman Hole East DDT 
Spitzer images. Five out of 10 submm sources are considered to 
have a secure radio detection (la RMS of 4.8 fiJy beam -1 , Ivi- 
son et al. 2002) and of these all are detected with IRAC and in the 
MIPS 24 pm data (3a depth of 120 /iJy, Egami et al. 2004). Using 
the same Spitzer data in the Lockman Hole East, Ivison et al. (2004) 
examined the Spitzer properties of MAMBO-selected galaxies and 
found counterparts for all 9 > 3cr MAMBO galaxies within the 
Spitzer field. Frayer et al. (2004) obtained targeted SCUBA pho- 
tometry of optically faint radio sources in the First Look Survey 
Verification (FLSV) field. Seven out of 28 sources are detected 
above 3a at 850 /im and all of these 7 are detected at 24 /j,m (al- 
though note that one source is only detected at 2.7a). In this paper 
we present the largest sample of submm galaxies with statistically 
secure Spitzer counterparts. 

In star forming galaxies, X-rays can come from star formation 
or weaker AGN activity (usually soft X-rays) or stronger AGN ac- 
tivity (usually hard X-rays). We have used the main and supplemen- 
tary catalogues from Alexander et al. (2003b) to look for any X-ray 
emission from the submm counterparts and attempt to determine 
whether the X-ray emission is due to AGN activity. We use the sim- 
ple criterion that any hard X-ray detected counterpart is considered 
an X-ray AGN. 11/21 (52 per cent) of the secure counterparts are 
X-ray detected in any band and 4/21 (19 per cent) are classified as 



X-ray AGNs. This AGN fraction is somewhat lower than that given 
in Alexander et al. (2005a, S^io per cent), although not incon- 
sistent, given the small number statistics. The AGN fraction could 
also be affected by variations in the X-ray and submm sensitivities 
across the field and/or differences in the derived AGN fraction from 
a purely submm- selected sample to that extrapolated from spectro- 
scopically identified submm galaxies. Both results are consistent 
with the majority of submm galaxies having their IR luminosity 
predominantly powered by star formation and not AGN activity. 
Although not as rigorous as the classification used in Alexander 
et al. (2003a) and Alexander et al. (2005b), in practice all of the 
AGNs in those studies were detected in the hard band. Through- 
out the figures in this paper, we denote submm counterparts which 
are classified as X-ray AGN (i.e. hard X-ray detected) with a cross 
symbol. We discuss the AGN signatures more in Section|5| 

4.1 Redshifts 

Table lA2l also provides redshift estimates for each source. When 
available we list the spectroscopic redshift, otherwise we list the 
photometric redshift (see Section lX2t . The redshift distribution of 
our sample is consistent with our previous redshift distribution pre- 
sented in Paper III and also with the sample of 73 radio-detected 
submm galaxies in Chapman et al. (2005). Note that 7/21 of our 
secure counterparts are also in the Chapman et al. 2005 sample 
and have spectroscopic redshifts. The median redshift for our se- 
cure counterparts with optical spectroscopic or photometric red- 
shifts (21/35 sources) is 2.0 (interquartile range 1.3-2.7), while the 
median redshift for all counterparts, including those with only IR 
photometric redshifts (33/35 sources), is 2.2 (interquartile range 
1.4-2.6). The large fraction of submm sources with reliable coun- 
terparts in our sample, which contains sources covering Ssso = 2- 
20mJy, seems to imply that there is not a large population of 
850 fim selected galaxies at z > 4. Even if we conservatively as- 
sume that one of the secure counterparts and two of the tentative 
counterpart are incorrect identifications and that they, and the two 
sources with multiple counterparts, are really at high redshift, this 
still constrains the fraction of submm sources at z > 4 to < 14 
per cent (< 5/35). There is evidence which suggests that galaxy 
samples selected at longer wavelengths (~ 1 mm) have a higher 
fraction of sources at z > 4 (Eales et al. 2003). 

4.2 New radio detections 

Using the published Richards (2000) 5a 1.4 GHz catalogue, we 
found radio counterparts for 37 per cent of our sample (see Paper 
II and III). We find that all of the VLA 1.4 GHz radio counterparts 
presented in Paper II and Paper III are confirmed in the new radio 
reduction and at similar flux levels. In addition, we find 18 new ra- 
dio sources above 3a in the new radio reduction within the search 
radius of a submm source, all of which are detected with IRAC 
and/or MIPS. The separation between the radio and IRAC positions 
of the submm counterparts is < 0.5 arcseconds. All of our secure 
counterparts are radio-detected and, including the tentative identifi- 
cations, 74 per cent (26/35) of submm sources in GOODS-N have 
a radio counterpart using the new VLA image. Note that 10 per cent 
of this improvement is simply due to lowering our signal-to-noise 
ratio threshold to 3a. Taking the 20 secure radio counterparts (ex- 
cluding GN14), 8/20 are only present in the new radio reduction. 
Radio flux densities for all of the sources detected above 3a are 
listed in Table lA2l and in the absence of a detection we list the 3a 
upper limit. 
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Submm-radio offset (arcsec) 

Figure 1. Submm-radio positional uncertainty in GOODS-N. The solid his- 
togram is the cummulatfve distribution of radial offsets between the submm 
and radio positions for the 21 secure radio counterparts in GOODS-N. The 
dashed curve is the expected radial offset distribution assuming a Gaussian 
with o r = a X (FWHM/SNR) submm , with a = 1.0. We have used 
the average submm SNR which has been corrected for flux deboosting. 
A KS test fails to find a significant difference between these two distri- 
butions. Ivison et al. (2005) found the same distribution for a robust sample 
of submm galaxies. The value of a found here is somewhat larger than what 
might be expected if centroiding errors dominate (e.g. Hughes et al. 1998). 
However, we expect astrometry shifts to be important in the heterogeneous 
GOODS-N super-map. Although the curves deviate at offsets greater than 
6 arcseconds, we do not feel justified in exluding them as reliable counter- 
parts, since they satisfy our criterion of having P < 0.05 and furthermore, 
3/4 of these counterparts are detected in other SCUBA photometry or map- 
ping data (Chapman et al. 2001, 2005). A uniform submm survey where all 
data were taken with the same chopping pattern, such as SHADES (Mortier 
et al. 2005), would be expected to provide a better test of the above relation 
(see Ivison et al. in preparation). 



Our confidence in these new radio detections is considerably 
bolstered by the fact that several of the radio sources correspond 
to the optical galaxies suggested as counterparts in Paper III, based 
only on optical/near-IR colours. The new radio counterparts that lie 
just above the detection threshold in this improved VLA reduction 
are consistent with our results of Paper III, which indicated that the 
radio-detected and radio-undetected submm sources are not drawn 
from very different populations. Extrapolating from this, we sug- 
gest that the remaining submm sources which are not currently de- 
tected in the radio are likely to lie just below the detection thresh- 
old. This is supported by the fact that several of them appear as 
~ 2a features in the radio image. 

Taking the 21 secure radio counterparts, we can investigate 
the distribution of offsets between the radio and submm position. 
Ivison et al. (2005) found that the distribution of offsets for robust 
submm sources detected at both 850 fim and 1.2 mm was approx- 
imately Gaussian with a ~ FWHM/SNR. In Fig.Q] we plot the 
cummulative distribution of radial offsets between the radio and 
submm positions for our 21 secure counterparts along with the 
expected Gaussian distribution. A Kolmogorov-Smirnov (KS) test 
does not find any significant differences between the two distribu- 
tions. 



4.3 Double radio sources 

When deep radio observations are present in a submm field, around 
10 per cent of the of the submm galaxies have two radio sources 
within the search radius (Ivison et al. 2002; Chapman et al. 2005). 
We find a consistent result in the GOODS-N field. With the new ra- 
dio reduction of the 1.4 GHz VLA data, we find that 4/35 submm 
sources have two radio counterparts in the search radius. Note 
that these additional radio sources were not present in the origi- 
nal Richards (2000) radio catalogue. The probability of having two 
radio sources of any brightness within 8 arcseconds of a submm 
source at random is less than 1 per cent, and therefore it is likely 
that the radio sources are associated with each other, and with the 
submm emission. Three of these double radio pairs (GN04, GN07 
and GN19) are very close (< 3 arcseconds) and two have an op- 
tical counterpart for only one of the radio sources, meaning that 
we are unable to get a redshift estimate of the second radio source 
(although all of the additional radio sources have a counterpart de- 
tected in IRAC and/or MIPS). As we discuss in the Appendix, if 
we assume both radio components are at the same redshift, then the 
separation between them is ~ 20 kpc, indicating that these could be 
interacting or multi-component systems. The IRAC data for these 
objects also support the notion that the radio sources are both at the 
same redshift and contributing to the total submm flux (see Section 
16.31 . Therefore we conclude that the double radio sources asso- 
ciated with GN04, GN07 and GN19 represent interacting systems, 
and hence we use the sum of the radio and mid-IR flux from the two 
components when determining the global multi-wavelength prop- 
erties of the submm source. 

For the other double radio source in our sample, GN17, both 
radio components have an infrared and optical counterpart and so 
we can estimate photometric redshifts, which we find to be different 
(•Zphot =1.7 and 1.2). In the next section, we describe how we fit 
the MIPS 24 fim flux and redshift to a suite of SED templates to 
determine the likely contributions to the 850 fim flux from each 
radio source. We found that, for GN17, the second radio/24 fim 
source is likely to be subdominant for the submm emission. 

GN38, which does not make our final submm catalogue due to 
flux boosting, has three radio sources within the search radius, two 
of which look to be associated. All of these double radio sources in 
GOODS-N are discussed individually in the Appendix, along with 
notes for the other submm sources. 



4.4 MIPS 24 fim detections 

The full width half maximum (FWHM) of the MIPS PSF is 5.7 
arcseconds, which is significantly better than the 14.7 arcsecond 
resolution obtained with SCUBA. In addition, detection of sources 
in the MIPS 24 image was carried out using the IRAC source 
positions as a prior, which facilitates matching between IRAC and 
MIPS images. There are a few cases where blending in the MIPS 
image is an issue, and we have dealt with these on a case by case 
basis to make sure we are getting the best estimate of the 24 fim flux 
at the position of the counterpart. As a result, all of the candidate 
24 fim counterparts within 8 arcseconds of the submm sources have 
a single IRAC counterpart. The measured FWHM of the IRAC PSF 
is < 2 arcseconds (Fazio et al. 2004), which is sufficient to identify 
unique optical counterparts to IRAC sources in most cases. 

Using no prior knowledge at any wavelengths, we have ex- 
plored submm counterpart identification using only the MIPS 
24 fim image. We find 10/21 secure counterparts have P < 0.05 
using only the MIPS 24 fim catalogue. These secure MIPS coun- 



© 0000 RAS, MNRAS 000, 000-000 



Spitzer's view of submm galaxies 1 



0.5 f 



0.4 



0.3 



0.2 



0.1 



Fraction with no MIPS in search radius: 
Submm positions -0.1 
Random positions -0.6 



I 



1.0 0.8 



0.6 0.4 
max S 24 /total S 24 



0.2 



0.0 



Figure 2. Multiple MIPS 24 fim sources within one SCUBA 850 fim beam. 
We have calculated the ratio of the total flux density of the brightest 24 fim 
source to the total flux of all 24 fim sources within a SCUBA 850 fim beam 
for our submm sources (shaded region) and for 1000 random positions in 
GOODS-N (dark line, slightly offset for clarity). These distributions have 
been normalized by the total number of sources in each sample. Note that 
a larger fraction of the random positions have no nearby 24 fim source. Of 
those positions which have at least one 24 fim source in the search radius, 
both distributions are consistent with being dominated ( > 50 per cent of 
flux) by a single 24 fim source most of the time. 



terparts are all brighter than 70 fiJy (> 14a) at 24 fim. While this 
is only 30 per cent of our original sample, it demonstrates that rel- 
atively shallow MIPS observations (5a ~ 70 /iJy) can be useful for 
identifying secure submm counterparts in the absence of radio data. 
Deep radio observations may be difficult to obtain in some fields, 
due to their location on the sky or the presence of bright quasars 
in the field which drive the dynamic range, and MIPS observations 
can be used as an alternative. In situations where we are fortunate 
enough to have deep radio and infrared observations, the high over- 
lap between 24 fim and submm sources provides an independent 
confirmation of the robustness of the radio counterparts. 

Roughly half (17/35) of the submm sources in GOODS-N 
have more than one > 3<r 24 fim source within the search radius. 
SCUBA observations are limited by confusion noise; there are 
many faint submm sources ( < 1 mJy) which get blended in the 
large JCMT beam. Hence it is possible that when there are mul- 
tiple 24 fim sources, they are all contributing to the submm flux. 
The important question is whether more than one of the mid-IR 
sources is making a significant contribution to the submm flux? We 
investigate this in two different ways. 

In the first test, we assume that the brightest 24 fim source 
within the SCUBA 850 fim beam will be associated with the 
submm source. While this may not always be the case, it is sta- 
tistically the most likely, given Poisson statistics. Note that we are 
not saying that we require the brightest 24 fim source to be the cor- 
rect submm counterpart (however, in practice this is usually true). 
In Fig. [2] we plot the distribution of the ratio of the total flux den- 
sity of the brightest 24 fim source to the total flux of all 24 
sources within a SCUBA beam (see also Sajina et al. 2006). This 
ratio provides a measure of the importance, in terms of 24 fim flux, 



of any additional sources in the search radius. A ratio close to 
one indicates that the submm emission is dominated by this sin- 
gle 24 fim source (under the assumption that 24 fim and 850 fim 
fluxes of sources are related). The results for our submm sources in 
GOODS-N are shown as the shaded region in Fig.|2|and, for com- 
parison, the dark curve is the distribution for 1000 random positions 
in GOODS-N. Roughly 60 per cent of the random positions do not 
have any 24 fim sources within the search radius, whereas fewer 
than 10 per cent of the submm positions have no 24 fim sources 
within the search radius. Of those positions which have at least one 
24 fim source in the search radius, both distributions are consistent 
with being dominated ( > 50 per cent of flux) by a single 24 fim 
source the majority of the time. Only about 10 per cent of the 
submm sources have more than one 24 fim galaxy with the addi- 
tional source(s) contributing up to ~ 50 per cent of the total MIPS 
flux in the SCUBA beam. Whether or not this additional source (or 
sources) contributes significantly to the submm flux depends not 
only on the MIPS flux, but also the SED and redshift, which we 
investigate next. 

In order to estimate how much the multiple 24 fim sources 
within a SCUBA beam are contributing to the total submm flux, we 
have fit the 24 fim flux of both (or three in some cases) sources at 
their redshifts (spectroscopic or photometric, see Section l4~l1 using 
a range of SEDs to determine the range of 850 fim flux expected for 
each source. Templates used include Chary & Elbaz (2001, here- 
after CEOl) and Dale & Helou (2002) star forming galaxy tem- 
plates, covering a wide range of luminosities, as well as an AGN 
template, Mrk231 (mid-IR spectrum from Rigopoulou et al. 1999, 
spliced with a fit to the Infrared Astronomical Satellite, IRAS, pho- 
tometry). We then see what combinations of sources produce the 
observed 850 fim flux of each submm system. Note that in this pro- 
cedure we are not assuming that the brightest 24 fim sources within 
the beam is the submm counterpart, since the fit will depend on 
the 24 fim flux and the redshift. Of course this procedure can only 
be a general indication, since we have not yet measured the entire 
SED shape of submm galaxies. However, we have been quite lib- 
eral in using a fairly wide range of templates, and it would take 
quite unusual SEDs to significantly change the results. Our conclu- 
sion is that, for the majority of the cases, one of the 24 fim sources 
is the dominant submm emitter (i.e. more than half of the submm 
flux comes from it). In most cases, the second 24 fim source in the 
SCUBA beam is usually faint and at low redshift, and therefore not 
capable (assuming the wide range of SED templates) of producing 
the observed submm flux. See the Appendix for detailed comments 
on individual sources. 



5 SED TEMPLATES 

The mid-IR portion of the SED of galaxies is sensitive to thermal 
dust emission, AGN power-law emission (Clavel et al. 2000) and 
polyaromatic hydrocarbon (PAH) emission features (Hudgins & 
Allamandola 2004). The MIPS 24 fim filter is narrow enough that 
strong spectral features have a significant effect on the photome- 
try. In particular, the 9.7 fim silicate absorption feature has been 
found to be very strong in IR-luminous galaxies (Dudley & Wynn- 
Williams 1997; Spoon et al. 2002, 2004a; Armus et al. 2004). The 
strength of this, and other, mid-IR spectral features in different 
galaxy types is still not completely understood, but deep IRS spec- 
troscopy of galaxies at low and high redshift is beginning to provide 
a more complete picture (Armus et al. 2004; Yan et al. 2005). 

Most models of the entire IR SED prior to Spitzer did not 
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properly characterize the mid-IR spectrum of galaxies, in particu- 
lar the effects on the broadband fluxes of varying PAH emission 
features and strong silicon absorption. CE01 provide SED tem- 
plates for galaxies as a function of infrared luminosity which re- 
produced the IR and submm observations of nearby galaxies at that 
time. More recently, Draine (2003 and references therein) has cre- 
ated a carbonaceous-silicate grain model which includes absorption 
features, such as the 9.7 fim silicate feature. This model provides 
wavelength-dependent extinction, which reproduces the observed 
obcuration of starlight and infrared emission. In order to account 
for possible effects on the mid-IR photometry from strong mid- 
IR spectral features, we have created a combined model which in- 
cludes the CE01 SED templates and additional extinction using the 
Draine (2003) grain models and absorption cross-sections, allow- 
ing for the amplitude of the additional extinction to vary. Hereafter 
we refer to these models as modified CE01 templates. 

At submm wavelengths, the SED is often simplified to a grey- 
body, whose shape is characterized by the temperature and /3 (Blain 
et al. 2002), which are related to the wavelength of the far-IR peak 
and the slope of the SED, respectively. The dust model for CE01 
was created from four components, at 18 K, 40 K, 300 K plus the 
PAH emission, and it is the relative contribution of each which 
varies as a function of luminosity. For our purposes, a cool CE01 
model has a higher contribution from the lower temperature com- 
ponents than the warm CE01 model, and therefore peaks at longer 
wavelengths. A rough estimate of the dominant temperature of the 
CE01 templates can be determined by considering the wavelength 
of the far-IR peak. In CE01, the temperature increases slowly with 
luminosity, due to the local observed luminosity-temperature rela- 
tion (see CE01 and references therein). However, at high redshift, 
little is known about the temperatures of IR-luminous galaxies, and 
so we have fit for luminosity (amplitude of SED) and temperature 
(shape of SED) separately in this paper. 

The radio portion of the CE01 models is estimated by assum- 
ing the local radio-IR correlation. Therefore we cannot say any- 
thing about whether or not the correlation holds for our sample, 
when we are using these models to fit for the IR luminosity. 

Arp220 is the best studied local ULIRG (Gonzalez- Alfonso et 
al. 2004; Spoon et al. 2004b) and is often compared to high redshift 
submm galaxies (e.g. Chapman et al. 2005). However, it is not rep- 
resentative of the entire sample of local ULIRGs, as its IR SED is 
cool compared to most IRAS-bright ULIRGs (Farrah et al. 2003), 
and has one of the most extreme mid- to far-IR slopes. We compare 
our results to the observed SED of Arp220 interpolated from pho- 
tometry listed in the NASA/IPAC Extragalactic Database. As we 
show in Section^ we find evidence that many high redshift submm 
sources have dust temperautres as cool as or cooler than Arp220. 
For completeness, we also compare to the observed SED of the lu- 
minous AGN Mrk231, where we have spliced the mid-IR spectra 
from Rigopoulou et al. (1999) with a fit to the IRAS photometry. 



6 INFRARED PROPERTIES 
6.1 Probing dust emission 

In high redshift submm galaxies, the 850 /im window is sensitive to 
cool dust close to the peak of the far-IR dust emission. On the other 
side of the far-IR peak, the mid-IR regime is also sensitive to ther- 
mal dust emission. However, as discussed in Section [5] there are 
possible additional contributions, including AGN power-law emis- 
sion and PAH emission features. As one moves to higher redshift, 
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Figure 3. The Sza/Sbbo ra ti° as a function of redshift for submm sources 
in GOODS-N. The open circles show the average values of the submm 
sources in five redshift bins, with roughly equal numbers of sources in each. 
The dotted is a cool SED template from CE01 models, the solid curve is 
the same CE01 template with additional additional extinction from Draine 
(2003), the dashed curve is the observed SED of Arp220 and the dash-dot 
curve is Mrk231, an IR-luminous AGN. The shaded region represents the 
redshifts where the 9.7 (im silicate feature passes through the 24 fim pass- 
band. The generally lower Sai/Sgso ratios for the submm sources suggest 
that they have higher levels of extinction, especially in the mid-IR, than the 
local ULIRGs which had been used to construct the CE01 models. 



the 850 fim filter climbs the smooth grey-body dust peak, while 
the 24 /im filter passes through various spectral features and also 
changes depending on the intensity of the AGN emission. For this 
reason, it would be surprising if there was a tight correlation be- 
tween S24 and 5850 as a function of redshift, although there should 
be a general trend. 

In Fig. [3] the SW/Ssso ratio is plotted as a function of red- 
shift for the sample of secure submm counterparts in GOODS-N. 
As listed in Table lA2l 17/21 secure counterparts have reliable red- 
shifts, with more than half of these being spectroscopic. Due to 
the many factors influencing the 24 fim photometry for individual 
sources (as discussed in Section[5}, we also plot the average ratio in 
4 redshift bins as the open circles. For comparison, we overplot sev- 
eral galaxy SED models: Arp220; Mrk231; a cool CE01 template; 
and the same CE01 template with additional extinction (see Section 
[5}. The shaded region in Fig.|3|denotes the area where the 9.7 fim 
silicate feature passes through the 24 fim passband. The cool CE01 
model (dotted curve) is a poor fit to the data and overpredicts the 
5 , 24/5 , 85o ratio at all redshifts. The solid curve in Fig.|5]is a modi- 
fied CE01 template (see Section|5J, where the additional extinction 
is of the order Av ~ 8 magnitudes (in Fig. [6] we will show that 
the total extinction (in the mid-IR) in this modified CE01 model is 
roughly equal to the extinction observed in Arp220 and therefore 
this model is entirely reasonable). Overall, the observed SW/Ssso 
ratio of the submm sources as a function of redshift appears consis- 
tent with the modified CE01 model. It also agrees reasonably well 
with the ratios expected for a redshifted Arp220 SED, particularily 
at higher redshifts. Submm sources classified as having an AGN 
due to the presence of hard X-rays do not have substantially higher 
•S'24/5'850 ratios, as one might expect if the mid-IR emission was 
dominated by the AGN. Overall, Fig. [3] suggests that the submm 
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sources have higher levels of extinction, especially in the mid-IR, 
than the local ULIRGs which had been used to constrain the CE01 
models. In particular, the silicate absorption feature at 9.7 /im may 
be attenuating the 24 /im flux. Based on 16/24 colours, Kasli- 
wal et al. (2005) predict that roughly half of ULIRGs at z ~ 1-2 
might be missed at 24 /im due to silicate absorption. As we dis- 
cuss in Section we require a cooler template when fitting the 
Ss5q/Siaghz ratio as a function of redshift. However, decreasing 
the temperature of the template is not enough to account for the 
low S24 1 Ssso ratios seen in Fig. [5] and therefore we require the 
additional extinction in the mid-IR. 

Another class of star forming galaxy at similar redshifts to 
the submm galaxies which has been well studied in the IR are the 
BzK galaxies (Daddi et al. 2004). These galaxies are selected on the 
basis of their stellar SEDs (through K-b&nA fluxes and optical-IR 
colours), rather than based on longer-wavelength (mid-IR, submm, 
radio) properties that trace dust heated by star formation, and there- 
fore they are subject to different selection effects than submm 
galaxies. While submm galaxies at 1.4 < z < 2.5 with a mea- 
surable K magnitude generally do satisfy the BzK criterion, the 
majority of submm counterparts at these redshifts have K > 21.9, 
and therefore there is very little overlap between galaxies selected 
at 850 /im and bright BzK galaxies. Only one (out of 21) of our se- 
cure submm counterparts qualifies as a K < 21.9 BzK galaxy but 
this increases to four if we also consider consider K > 21.9 BzK 
galaxies. Daddi et al. (2005) found that 82 per cent of the GOODS- 
N BzK galaxies with K < 21.9 (AB magnitudes) are detected at 
24 /im, with an average flux density of about 125 /ily. At 2 ~ 2, 
they have typical S24 / Ssso ratios of 0. 1 6 (Daddi et al. 2005), which 
is larger (by a factor of 4) than that of most of the submm sources 
at z ~ 2. This might indicate that their SEDs peak at shorter wave- 
lengths. However, a difference in AGN contribution and PAH fea- 
tures in the IR luminosity could also cause this effect. BzK galax- 
ies (K < 21.9) are ~5 times less luminous than typical submm 
galaxies (see Table A2 and Chapman et al. 2005). The higher lumi- 
nosities of 850 /im selected galaxies (and lower surface densities) 
means that, at best, submm sources represent a small fraction of 
the most extreme BzK galaxies (Daddi et al. 2005; Dannerbauer et 
al. 2006). 

6.2 Probing stellar emission 

At the redshifts of these SCUBA sources, the IRAC channels are 
sampling the rest-frame near-IR, which is sensitive to direct emis- 
sion from stars. If these galaxies are not dominated by an AGN at 
these wavelengths then the IRAC colours should show an indication 
of the 1.6 /tm local peak of the f v stellar SED. In Fig.|4] we plot 
an IRAC-based colour-magnitide diagram, specifically 5*5. 8/S3. 6 
as a function of S5.8 for the submm sources and the GOODS-N 
field galaxies. The redshift of the submm source is indicated by 
the symbol. We see that z ~ 1.5 provides a clear transition in 
this IRAC colour for the submm sources, due to the filters moving 
over the stellar peak. The sources which lack redshift information 
are all consistent with the colours of the higher redshift submm 
sources. These sources are all faint in the optical, which indicates 
that they are either at high redshift or have a very high dust extinc- 
tion (or both). For a given colour, the submm galaxies are brighter 
at 5.8 /im, which is sensitive to stellar mass, indicating that they 
are on average more massive (Borys et al. 2005). The presence of 
an AGN could complicate matters, and so we need to check for 
any bias. Sources which are detected in the hard-band X-ray image 
are indicated in Fig.|4]by a cross, and there seems to be no strong 
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Figure 4. IRAC-based colour-magnitude diagram for GOODS-N. Larger 
symbols are our submm sources, while the small black dots are field galax- 
ies which are detected in both IRAC bands. Filled circles and squares denote 
submm sources with spectroscopic and photometric redshifts, respectively. 
Crosses are sources which contain an AGN, as indicated by the presence 
of hard X-rays. The black solid symbols are at low redshift (2 < 1.5), the 
grey solid symbols are at high redshift (2 > 1.5) and the black open sym- 
bols have unknown redshift. Flux uncertainties are shown as error bars, al- 
though they are smaller than the symbol size at brighter fluxes. The submm 
galaxies are bright at 5.8 /tm for a given colour compared to the field galax- 
ies. Since 5.8 /tm samples near the peak of the stellar SED in these sys- 
tems, this could mean that the submm galaxies are on average more massive 
(e.g. Borys et al. 2005), although the presence of an AGN will also affect 
the rest-frame near-IR flux. This IRAC colour clearly separates the submm 
sources into low and high redshift sub-samples, due to the fact that these 
wavelengths are sampling the peak in the stellar SED. This figure suggests 
that the sources with unknown redshift are at 2 > 1.5, due to their steep 
S5.8/S3.6 colour. The dashed lines indicate the region where S5. s > 6 
and S5.8/S3.6 > 1-3, within which 11/11 of the radio-detected submm 
counterparts are confirmed to be at 2 > 1.5. This IRAC colour-magnitude 
cut can be used to help identify submm counterparts, since there is a low 
probability of landing in this region of the diagram at random. 



colour distinction between the sources with and without a hard X- 
ray detection, although the former are generally brighter at 5.8 /im. 
BzK galaxies (with K < 21.9) are also bright at 5.8 /im, indicating 
that they are massive. However, their S5. s/Ss.e colour is slightly 
bluer than that of the submm galaxies of similar redshifts. This may 
indicate that, while this IRAC colour is primarily sensitive to red- 
shift, the Ss.s/S^s colour of submm galaxies may also be affected 
by an AGN or dust obcuration. 



6.3 Spitzer photometric redshifts 

Expanding on the separation of IRAC colour with redshift seen in 
Fig.|4| we have attemped to estimate IR photmetric redshifts for this 
sample using only the IRAC and MIPS photometry. The idea is to 
solve for the redshift as a linear combination of the logarithm of the 
flux density at the mid-IR wavelengths (see Equation|2j to provide a 
completely model independent estimate of the redshift (Connolly et 
al. 1995; Sajina 2006). We use the secure submm counterparts with 
known optical spectroscopic redshifts to solve for the coefficients 
and then apply the formula to sources with unknown redshifts. This 
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Figure 5. IR photometric redshift accuracy. We plot the (model indepen- 
dent) IR photometric redshift, as determined only from IRAC and MIPS 
photometry, versus the known optical spectroscopic redshift for 10 secure 
submm counterparts in GOODS-N. Crosses indicate counterparts which are 
classified as X-ray AGN due to the presence of hard X-rays. Spitzer pho- 
tometry can be used to constrain the redshift for sources which are faint or 
undetected in the optical. While the IR photometric redshifts are not as ac- 
curate as the optical photometric redshifts (see Pope et al. 2005), they are 
good to Az < 0.4 for 8/10 sources. 



is particularily useful for the submm counterparts which are faint 
or undetected in the optical and thus lack an optical photometric 
redshift estimate. The redshift is calculated using the following for- 
mula: 

ZIR = a + b ■ log(53.6) + c ■ log(5 5 .8) 

+ d ■ log(&.o) + e • log(&4). (2) 

Using the 10 secure submm counterparts with spectroscopic 
redshifts, we fit for the five coefficients which give the lowest RMS 
when compared to the spectroscopic redshifts. Fig.|5|shows the ac- 
curacy of the IR photometric redshift for these 10 submm sources. 
The RMS dispersion of the redshift errors for these submm galaxies 
is cr(Az/(l + z)) = 0.07, using values of 3.3, -2.5, 4.6,-1.1 and 
-1.4 for a, b, c, d and e, respectively. We have performed this fit for 
various combinations of these flux densities, (including functions 
without the logarithm), and found this combination to give the low- 
est RMS for this sub-sample of galaxies (although it is by no means 
unique and different coefficients could also give low RMS). Using 
the 4.5 pirn flux density instead of (or as well as) the 24 flux 
density resulted in an RMS which was higher, and a coefficient for 
the 5*4.5 term which is close to zero. This is not surprising because 
4.5 fj,m samples the peak of the stellar light distribution at z ~ 1.8, 
and since we are not fitting for any sources near this redshift, the 
inclusion of this parameter does not improve the redshift estimate. 

For the submm counterparts which lack an optical photomet- 
ric redshift estimate we list the zm. in Table lA2l in brackets. While 
these IR-derived redshifts are not as accurate as optical photomet- 
ric redshifts for submm galaxies (Paper III), they do provide an 
accuracy of Az < 0.4 for 8/10 secure counterparts with spectro- 
scopic redshifts. This is not accurate enough for exact determina- 
tion of physical parameters, but it can be useful for separating the 



high redshift and low redshift sub-populations of submm galax- 
ies, particularly for sources where the optical counterpart is invis- 
ible or ambiguous. This technique only works for particular sub- 
samples of galaxies (i.e. Equation[2|does not apply to other galaxy 
samples), since we are essentially assuming that the SEDs have a 
similar shape at these wavelengths (i.e. fitting for the coefficients 
using all spectrscopically confirmed IRAC sources in GOODS-N 
would give less accurate IR photometric redshifts for the submm 
sources). However, this technique will likely improve with larger 
submm samples with known redshifts. As we learn more about the 
range of SEDs for submm galaxies, then photometric redshift esti- 
mators using the entire wavelength range should be able to acheive 
much better results. 



7 SPECTRAL ENERGY DISTRIBUTION 

The addition of the IRAC and MIPS photometry to the multi- 
wavelength dataset for submm galaxies provides powerful con- 
straints on the shape of the SED, the nature of the power source, 
and the source redshift. In Fig. |6| we plot composite rest-frame 
SEDs. All photometry points from the near-IR through to the ra- 
dio are plotted, with each colour indicating a different passband. 
Note that we have only included the submm sources with spectro- 
scopic redshifts, since we would like to isolate the scatter due to 
different SEDs and/or luminosities; if we also include the sources 
which have photometric redshifts, it does not qualitatively change 
these figures, but the scatter is increased. By minimizing the scat- 
ter (RMS of the logarithm of the 24 fim, 850 /xm and 1 .4 GHz flux 
densities relative to the model), we find the best-fit modified CE01 
template (see Section|5J. This model (solid curve in Fig.[6j has an 
IR luminosity of 3 x 10 12 Lq. The additional extinction suppied 
by the Draine (2003) models amounts to Av ~ 8 magnitudes. In 
the top panel, we plot the rest-frame composite SED in physical 
units. In the middle, and bottom, panels, the SED is normalized 
to the submm and radio flux densities, respectively, to match the 
modified CE01 template. The three panels of Fig. [6] can be used 
to assess how well this model fits the submm sources using these 
different normalization approaches, as we discuss below. 

The main conclusion to be drawn from this plot is that high 
redshift submm galaxies appear to have different SED shapes from 
those observed locally for galaxies of similar luminosities. A major 
contribution to this difference could be the temperature of the dust, 
but it is important to realize that this just means an SED peaking at 
longer wavelengths. The best-fit submm SED (solid curve in Fig.|6j 
peaks at ~ 100 /im (T ~ 30 K), while a typical local ULIRG tem- 
plate from CE01 (dash-dot curve in top panel of Fig.[6) of the same 
luminosity peaks at ~ 85 /im (T ~ 34 K). The same point is illus- 
trated in Fig.0 where we plot the Sg50 /Si.4GHz ratio as a function 
of redshift. Again, the submm counterparts are inconsistent with the 
warm local ULIRG template and fit better to the cooler template. 

Cooler temperatures for the submm galaxies is consistent with 
predictions from hierarchical galaxy formation models (Kaviani, 
Haehnelt, & Kauffmann 2003). Similar results are also seen in 
Blain et al. (2004) and Chapman et al. (2005) for a large sample 
of 73 bright radio-detected submm galaxies. Using the redshift, ra- 
dio and submm flux, Chapman et al. (2005) found that the median 
dust temperature of their sample is lower than that estimated for 
local IRAS 60 /im galaxies of the same luminosity. The SCUBA 
Local Universe Galaxy Survey (SLUGS, Dunne et al. 2000; Dunne 
& Eales 2001) measured the submm fluxes of galaxies from the 
IRAS-bright galaxy sample. They find an average dust tempera- 



© 0000 RAS, MNRAS 000, 000-000 



Spitzer's view ofsubmm galaxies 1 1 



10 2 



N 

5 io 25 



g 10 i4 r 



3 



10 2 



10 2 




I I I 1 1 1 1 — I 1 1 1 I — I — I 1 1 1 m l — I 




1.0000 



0.1000 r 



T3 

§ 0.0100 



Z 0.0010 



0.0001 



100.00 r 



10.00 



1.00 



0.10 



0.01 

10 1 10 2 10 3 10 4 
Rest wavelength (|xm) 

Figure 6. Composite rest-frame SEDs for the GOODS-N submm sources. 
IRAC through radio photometry data are plotted, with each observed-frame 
wavelength shown with a different colour. In all panels, we only plot sources 
with spectroscopic redshifts since we are trying to isolate scatter due to lu- 
minosity and/or SED shape. The top panel shows the SED in physical units, 
while in the other two panels we have normalized the distribution to the 
submm and radio flux, respectively. In all panels, the solid curve is a mod- 
ified CE01 model for a cool starburst galaxy, where we have applied addi- 
tional extinction using the Draine (2003) models. This template has a total 
IR luminosity (8-1000 /ton) of 3 X 10 12 Lq and a temperature of 30 K. In 
the top panel, the dash-dot curve is a CE01 model with the same luminosity 
but a shorter peak wavelength. The dashed curve in the middle and bottom 
panels is the observed SED of Arp220. Error bars on the photometry are 
typically much smaller than the scatter among the sources, and have been 
omitted for clarity. 
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Figure 7. The Ssso / Sughz rat '° as a function of redshift for submm 
sources in GOODS-N. Filled circles and squares denote submm sources 
with spectroscopic and photometric redshifts, respectively. The crosses de- 
note sources which contain an AGN, as indicated by the presence of hard 
X-rays. The open circles are the average values of the submm sources in 
four redshift bins with roughly equal numbers of sources in each. The solid 
curve is a cool SED template from CE01 models, while the dashed curve is 
a warmer CE01 template and the dash-dot curve is Mrk231, an IR-luminous 
AGN. The average 5850 / Si.4GHz ra ti° f° r the submm sources suggest that 
they have SEDs which are cooler (i.e. peak at longer wavelengths), particu- 
larily those sources at lower redshifts. This may be the result of the selection 
at 850 fim. 



ture of 35.6 ± 4.9 K. However, their sample is less luminous than 
the blank sky submm galaxies, with IR luminosities in the range 
10 10 L Q to 3 x lO n i . Farrah et al. (2003) fit for the dust tem- 
peratures for a sample of local ULIRGs including faint IRAS galax- 
ies, and find a median temperature of 32 K for sources with an IR 
luminosity > 1O 12 L0. 

This result emphasizes the strong selection effects, both lo- 
cally and at high redshift. Our knowledge of local ULIRGs is 
not immune to selection effects, since it is dominated by results 
from IRAS which may preferentially select galaxies with warmer 
SEDs (peaking at shorter wavelengths). Similarity at high red- 
shift, 850 imaging preferentially selects galaxies which peak 
at longer wavelengths, and therefore comparing the two becomes 
problematic. Our average value of 30 K is slightly lower than the 
36 K estimated in Chapman et al. (2005) for brighter submm galax- 
ies, although this is probably due to the fact that our sample con- 
tains more low redshift, less luminous objects, where the 850 
selection effects are even stronger (see Fig.0. The submm selec- 
tion effects also mean that we may be missing a population of hot- 
ter ULIRGs which could contribute signitificantly to the star for- 
mation at high redshift (e.g. Chapman et al. 2004b). We require 
more objects selected at shorter submm wavelengths and longer IR 
wavelengths to test this hypothesis and determine if there are actual 
differences in the SED shape of ULIRGs at high and low redshift, 
or if we are just missing a large fraction of ULIRGs at high redshift. 

If high redshift star-forming galaxies really do have cooler 
temperatures then this may indicate that their emission is more 
extended than that of local ULIRGs, in which the majority of 
the IR emission comes from within the central kpc (Charman- 
daris et al. 2002). In support of this, Chapman et al. (2004a) 
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Table 2. RMS scatter in the SEDs of submm galaxies relative to the mod- 
ified CE01 model. Fig.lolshows the composite SEDs for spectroscopically 
identified submm galaxies in GOODS-N. In this table, we quantify the scat- 
ter (RMS of the logarithm of the flux densities) seen in the data points with 
respect to the best-fit model (cool CE01 model with additional extinction 
from the Draine 2003 models). The rows list the scatter (in dex) as a func- 
tion of observed wavelength. The columns list the scatter when the SED is 
in physical units, normalized to the submm, and radio, respectively (corre- 
sponding to the three panels in Fig.l6l. 
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found that 8/12 submm galaxies were spatially extended (on scales 
of ~ lOkpc) in their radio emission (see also Ivison et al. 2002 
and Muxlow et al. 2005), however much larger samples of radio- 
detected submm galaxies are needed in order to come to any defini- 
tive conclusions. 

Focussing now on the middle panel of Fig.|6| by normalizing 
to the submm flux, we are essentially assuming that Sssoum is a 
proxy for the far-IR luminosity. This makes sense, since the nega- 
tive K-correction almost cancels distance dimming over the typical 
z ~ 1-4 redshift range (see e.g. Blain et al. 2002). But if this is 
strictly true then the SED normalized at 850 /im will be indepen- 
dent of luminosity and should show variations in the SED shape 
only. 

In Table [2] we quantify the scatter (RMS of the logarithm of 
the flux density) at 24 /im, 850 /im and 1 .4 GHz for all three panels 
in Fig.|6| relative to the modified CE01 model. We have determined 
that the scatter due to photometric errors is negligible at this level, 
and therefore the scatter in the panels of Fig.|6|represents ranges in 
IR luminosity and SED shape. The scatter at 24 fim and 1 .4 GHz is 
lower for the submm normalized SED than for the physical SED, 
which implies that on average 850 /im is a good proxy for luminos- 
ity. In addition, the fact that there is still scatter after normalizing 
to either submm or radio flux (i.e. after removing the scatter due 
to luminosity) implies that the variations in SED shape are con- 
tributing significantly to the overall scatter seen in the top panel 
of Fig. [6] While we have not used the IRAC points to fit to the 
model SED, we notice that the intrinsic scatter of the IRAC points 
is lower when normalized to submm flux. Since IRAC is sensitive 
to stellar mass and 850 /im is sensitive to dust mass at the redshifts 
of these sources, this might indicate a correlation between stellar 
mass and the dust mass. In order to test this further requires indi- 
vidual optical-radio SEDs fits for a large sample of submm galaxies 
and we reserve this for a future paper. 



8 INFRARED LUMINOSITIES AND STAR FORMATION 
RATES 

8.1 Estimating LIR and SFR 

Given that the average submm galaxy seems to fit well by the mod- 
ified CE01 templates when both temperature and luminosity are al- 
lowed to vary, we have estimated the IR luminosity of each submm 
source by fitting its 24 /im, 850 fim and 1.4 GHz flux density to a 
suite of modified CE01 templates. In particular, we are interested in 



determining which of the three wavelengths is a better indicator of 
luminosity and thus star formation rate (SFR) in the submm galax- 
ies using current templates. We have fit for IR luminosity (LIR) us- 
ing the redshifted flux density at each wavelength separately (giv- 
ing us LIR24, LIR850 and LIR1.4), and using all three together (giv- 
ing us LIRall). We use the flux uncertainties as weights when fit- 
ting for the best template and checked to make sure that all data 
points (including the uncertainty) were consistent with the best- 
fit template. When we are fitting for LIR24, LIRgso or LIR1.4, we 
only have one data point, so we cannot constrain luminosity (nor- 
malization of the SED) and temperature (shape of the SED), and 
therefore we must assume that the CE01 templates accurately rep- 
resent the luminosity-temperature relationship and fit for just the 
luminosity. However, when we are fitting all three data points to 
get LIRall, we can allow both temperature and luminosity to vary 
independently. We consider LIRall to be the best luminosity esti- 
mator, since it uses the most information, has the widest range of 
templates, and fits best to all three data points. 

We define IR luminosity as the integral under the SED from 
8-1000 /im. We have fit for the IR luminosity of all submm 
sources which have a redshift estimate (either photometric or spec- 
troscopic) and the results are shown in Fig. [3] Using templates 
based on local ULIRGs with a fixed luminosity-temperature re- 
lation (CE01), the submm flux will typically overestimate the lu- 
minosity. This effect is most dramatic for lower luminosity submm 
objects which, as we will see in the next section, tend to lie at lower 
redshifts where the 850 /im window samples even further from the 
IR SED peak. However, when the temperature and luminosity are 
allowed to vary, 850 /im, 24 /im and 1.4 GHz flux density together 
give a good fit to the luminosity. Both 24 fim and radio flux ap- 
pear to be good indicators of the IR luminosity on their own, re- 
gardless of the 850 /im flux. This is because the rest-frame near- 
IR and radio are less sensitive to temperature variations. Therefore, 
submm galaxies appear to follow the local mid-IR radio correlation 
but not the local far-IR (200-400 /im) radio correlation, since they 
appear to have cooler temperatures than local ULIRGs (as discussd 
in Section^. Reddy et al. (2006) also find that the infrared lumi- 
nosity, as estimated from the submm flux alone is higher than that 
estimated using only 24 /im flux for a sample of 9 radio-detected 
submm galaxies with spectroscopic redshifts. 

Dale et al. (2005) predict a very large scatter in the IR luminos- 
ity as estimated at 24 /im (rest-frame 8 /im for the average submm 
source) in high redshift galaxies, based on observations of rela- 
tively local galaxies in the Spitzer Infrared Nearby Galaxies Sur- 
vey (SINGS, Kennicutt et al. 2003). These galaxies are much less 
luminous than submm- selected galaxies and are perhaps forming 
stars under different physical conditions than starbursting galaxies 
at high redshift. In our sample of 850 /im selected galaxies, we do 
not see the large scatter in the IR luminosity predicted by Dale et 
al. (2005), however it could be that we are preferentially only se- 
lecting the cool ULIRGs at high redshift. The mid-IR spectral fea- 
tures, including the 9.7 /im absorption feature, do not seem to have 
a large effect on the estimate of the luminosity. Also the presence 
of an X-ray counterpart does not appear to change the IR luminos- 
ity estimated by the 24 /im, 850 /im and 1.4 GHz flux, as might be 
expected if the AGN was strongly influencing the mid-IR spectra. 

Daddi et al. (2005) found that the mean 24 /im, 850 /im and 
1.4 GHz fluxes for the K < 20 BzK sample agree within a factor 
of 2 with those predicted by an unmodified CE01 template with IR 
luminosity of 1-2x10 12 Lq. While their redshift ranges overlap 
significantly, submm galaxies in our sample have average 24 /im, 
850 /im and 1.4 GHz fluxes which are a factor of 2, 7 and 2 times 
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Figure 8. Infrared luminosities and SFRs for submm sources detected at both 24 (im and 1.4 GHz. Luminosities are estimated by integrating from 8-1000 /jm 
under the best-fit modified CE01 template. The bottom axes show the luminosities estimated using only one wavelength and the redshift. The left-hand axis 
is the estimated IR luminosity using data at all three wavelengths. Filled circles and squares denote submm sources with spectroscopic and photometric 
redshifts, respectively. The crosses denote sources which contain an AGN, as indicated by the presence of hard X-rays. The top and right-hand axes show the 
corresponding SFR using the expression in Kennicutt (1998). When only one data point is available to constrain the SED then the IR luminosity is overpredicted 
by 850 fim, since the temperature is not allowed to vary and it is forced to fit to a warmer template. However, when the temperature and luminosity are allowed 
to vary, 850 fim, 24 (im and 1.4 GHz flux densities together give a good fit to the luminosity. 24 fim and 1.4 GHz data are not as sensitive to temperature 
and therefore they are able to predict an IR luminosity which is consistent with the result obtained when fitting all three and allowing both temperature and 
luminosity to vary. 



larger, respectively, than those for the BzK galaxies in Daddi et 
al. (2005). BzK galaxies seem to have similar SED shapes (temper- 
atures) to local ULIRGs, while the SEDs of submm galaxies peak 
at longer wavelengths. Submm galaxies appear to be quite atypi- 
cal (being both more luminous and cooler) when compared to local 
ULIRGs and high redshift BzKs. Vega et al. (2005) have recently 
suggested that enshrouded star forming galaxies undergo four dis- 
tinct phases characterized by different mid- and far-IR colours. Un- 
der their scheme, our results suggest that the submm galaxies, with 
cooler temperatures, may be an earlier phase in the star formation 
process and the BzK galaxies are a later, more evolved, phase (see 
also Dannerbauer et al. 2006). 

Individual values for the IR luminosity of the submm sources 
with spectroscopic redshifts are listed in Table A2. These are the 
luminosity esimates from the best-fit template using all three flux 
densities and the redshift. The mean and median IR luminosity is 
6.7 and 6.0 x10 12 Lq. These values are consistent with the lu- 
minosity values derived in Chapman et al. (2005) by fitting the 
Ss5q/Sia for 73 radio-detected submm galaxies to the Dale & 
Helou (2002) templates (although their definition of IR luminos- 
ity is slightly different, extending to 1100 fim). 

In Fig. [8] the top and right-hand axes give the corresponding 
value of the SFR, following the relationship between SFR and IR 
luminosity for starburst galaxies given in Kennicutt (1998): 

SFR (Mq yr -1 ) = 1.8 x 10" 10 L 8 _i oo , m (L ). (3) 

This relation assumes a Salpeter (1955) intial mass function and ap- 
plies to starbursts with ages less than 100 Myr. It also assumes little 
or no AGN contribution to the IR luminosity (see Section [5). The 
mean and median SFR for this sub-sample of submm sources in 
GOODS-N is 1200 and 1100 Mgyr" 1 , respectively. This is con- 
sistent with previous studies, telling us that submm galaxies are 
significant contributors to the global star formation at high redshift 
(e.g. Hughes et al. 1998; Lilly et al. 1999; Barger et al. 2000). How- 
ever, this is the first time that such values have been estimated using 
the mid-IR as well as the radio and submm regions of the SED. 



8.2 Evidence for evolution? 

Although the IR luminosity derived from the submm flux (left panel 
of Fig. |8j is higher than that derived from all three wavelengths, 
there is still a correlation between the LIRsso and LIRall- To 
investigate this, in Fig. [5] we plot the luminosity (LIRall) sepa- 
rately as a function of the 850 fim flux density and the redshift. 
Note that plots of the 3.6, 4.5, 5.8, 8.0, 24 fim and 1.4 GHz flux 
density as a function of 850 fim flux density do not show any sig- 
nificant trends and are dominated by the scatter. The correlations 
seen in Fig.|5|were already suggested in Paper III, where we found 
a lack of submm-faint objects at high redshift (see Fig. 3 of Pope et 
al. 2005). We concluded that there might be several factors affect- 
ing this trend, including the strong evolution of the number density 
of ULIRGs with redshift. Ivison et al. (2002) also pointed out the 
tendency for brighter submm sources to lie at higher redshifts and 
the recent results of Wang, Cowie & Barger (2005), which show 
that a large fraction of the 850 fim backgroud light comes from 
sources at low redshift (z — 0-1.5), are also consistent with this. 
The right panel of Fig. [5] which shows that the IR luminosity, as 
determined by the mid-IR, submm and radio flux (LIRall), is also 
correlated with redshift, strengthens the case that we are seeing the 
effects of evolution. At higher redshifts, the lack of lower lumi- 
nosity 850 fim selected sources is surprising, since the increase in 
volume implies that we should see roughly twice as many as those 
seen at lower redshifts. While it is contributing, the effect of volume 
cannot completely account for the trends seen in Fig. [5] 

The fact that the submm galaxies are characterized by cooler 
SEDs than local ULIRGs of the same luminosity implies that the 
typical SED temperature of ULIRGs may change with redshift, 
which could contribute to the trend in the left panel of Fig.|5| Lower 
luminosity submm galaxies at z > 2 would have to be character- 
ized by an SED that peaks at even longer wavelengths in order to be 
detected at 850 fim. It is worth noting that the sample of 73 radio- 
detected submm galaxies in Chapman et al. (2005) also shows a 
lack of low luminosity objects (LIR < 2 x 10 12 Lq) at high redshift 
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Figure 9. IR luminosity as a function of submm flux density and redshift. The quantity on the j/-axis is the 8-1000 fim IR luminosity derived from fitting 
the 24 fim, 850 fim and 1.4 GHz flux densities to the CE01 templates. Filled circles and squares denote submm sources with spectroscopic and photometric 
redshifts, respectively. Note that only submm sources with secure counterparts are plotted in this figure. The crosses denote sources which contain an AGN, 
as indicated by the presence of hard X-rays. In the left panel, we plot lines of constant LIR (1O 12 L0)/Sg5o = 0.34, 0.86, and 1.3, which correspond to the 
lower quartile, median and upper quartile for this sub-sample of submm sources. While there is a linear correlation in these data, the scatter (roughly two 
orders of magnitude) indicates that the submm sources are characterized by a wide range of SED shapes (since the negative K-correction largely removes the 
effects of redshift). The right panel shows that the IR luminosity of submm counterparts is a strong function of the 850 ^m flux. A plot of redshift as a function 
of 850 fim flux for this sample is presented in Pope et al. (2005) and also shows a correlation. These plots are difficult to understand without invoking the 
evolution of ULIRGs with redshift. 



(z > 1.5). To fully test the conclusions hinted at by our sample and 
that of Chapman et al. (2005) will require thousands of submm se- 
lected galaxies. Surveys with SCUBA-2 will produce these large 
samples and, along with extensive multi-wavelength data, will al- 
low for a detailed study of the luminosity function of 850 /im se- 
lected galaxies. 



9 AGN CONTRIBUTION 

In the previous section, we assumed that IR luminosity and SFR 
are directly related. But are the very large IR luminosities of 
submm galaxies powered by starbursts, active nuclei, or both? 
Many submm galaxies show evidence of an AGN, as determined ei- 
ther through X-ray observations (Alexander et al. 2003b, 2005b), or 
with optical spectra (Ivison et al. 1998; Swinbank et al. 2004; Chap- 
man et al. 2005). However, the most important question from the 
point of view of understanding the galaxy properties is not whether 
there are detectable AGN signatures, but whether or not the AGN is 
a significant contributor to the bolometric luminosity of the galaxy. 
Although, this question is difficult to answer definitively, data in 
the mid-IR can help determine what powers submm galaxies. 

With deep IRS spectra, Lutz et al. (2005) have determined the 
relative AGN contribution for two well-studied submm galaxies. 
One galaxy in the Lutz et al. (2005) sample shows that the AGN 
is a minor contributor to the luminosity, while the other seems to 
be powered equally by a starburst and an AGN. This small sample 
was limited to spectroscopically identified submm galaxies which 
are bright at 24 /im, and therefore it seems reasonable to expect 
that the whole submm population will show an even wider range of 
AGN properties. Both Egami et al. (2004) and Ivison et al. (2004) 



Table 3. Candidate AGN-dominated submm sources in GOODS-N from 
our secure counterpart list, either hard X-ray detected (Alexander et 
al. 2005b), or with a power-law shape in their mid-IR SEDs. The second 
column indicates if the source is detected in any X-ray band, while the third 
columm indicates if it is hard X-ray detected. For the final column we as- 
sume that the mid-IR SED has the form L„ oc u~ a and fit the IRAC data 
points to this model for each submm source. If the fit is good (low \ 2 ) then 
we list the best-fit value for a in the last column. Based only on these two 
AGN indicators, hard X-rays and mid-IR SED, we find that there is only 
one submm source in this list, GN04, which is likely to be dominated by 
the AGN in the mid-IR. 
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have used the mid-IR photometry as a diagnostic tool to put limits 
on the contribution made to the IR luminosity by an AGN and agree 
that this contribution is < 25 per cent for the SCUBA/MAMBO 
sources. Alexander et al. (2005b) put constraints on the AGN con- 
tribution to the energetics of a larger sample of radio-detected 
submm galaxies on the basis of their absorption-corrected X-ray lu- 
minosities. They found that, on average, the AGN contribution was 
likely to be negligible (~ 10 per cent) unless submm galaxies have 
a substantially larger dust-covering factor than optically selected 
quasars. The data from mid-IR and X-ray observations are comple- 
mentary and both are required to place the tightest constraints on 
the power source of submm galaxies. 

Most of our sample is quite faint at 24 pm, and so IRS spec- 
troscopy is not yet available. As shown in Fig. [3] the AGN Mrk231 
has a higher Sm/Ssso ratio, due to additional contribution to the 
mid-IR flux from the warm dust which is heated by the AGN. The 
low 5 I 24/S , 85o ratios for the submm galaxies argue that an AGN 
does not dominate the bulk of the IR emission in the submm sys- 
tems. We can use the optical through mid-IR photometry to look for 
a significant AGN power-law component. In addition, we can use 
hard X-rays to indicate the presence of an X-ray AGN, although 
a lack of hard X-rays does not rule out the possibility of an AGN 
being present. Note that it is important to look at both indicators 
of AGN presence, since locally we see powerful AGN which look 
like a starburst in the infrared while their X-ray properties reveal 
the AGN (e.g. NGC6240, Vignati et al. 1999; Lutz et al. 2003). On 
the other hand, we also see cases where the mid-IR SED does a 
better job of identifying the AGN and the X-rays are inconclusive 
(e.g. NGC1068, Jaffe et al. 2004). 

In Table [5] we summarize the candidate AGN-dominated 
submm sources in our sample. All sources listed here are either X- 
ray detected or have a power-law spectrum in the optical/IR regime. 
We find that 3/21 secure submm counterparts (14 per cent) emit 
hard X-rays, indicative of an AGN. For 4/21 sources (19 per cent), 
the colours in the IRAC bands can be well fit by a single power law 
L„ oc v~ a , whose spectral index a we report in Table|5| However, 
even if the mid-IR SED appears to follow a power law, this does not 
guarantee that the IR luminosity is dominated by the AGN, since 
there could be non-AGN components contributing to the mid-IR 
SED. We find that, using both X-ray and mid-IR data, there is only 
1/21 (5 per cent) source in our sample which has the potential for 
having its IR luminosity dominated by an AGN: GN04. Note that 
this is only an indication of the fraction of AGN-dominated sources, 
since we are not probing the far-IR luminosity directly. The com- 
bination of mid-IR spectroscopy and X-ray information will im- 
prove the census of AGN activity in these sources and will pro- 
vide the most definitive evidence on what is powering the energet- 
ics in submm galaxies, as well as helping to unravel whether there 
is an evolutionary connection between luminous submm galaxies 
and quasars. 



10 CONCLUSIONS 

We present multi-wavelength identification of 33 submm galaxies 
in GOODS-N using the 5p;fzer/GOODS Legacy data and a new re- 
duction of the VLA 1 .4 GHz radio data. 2 1 of these identifications 
have a low probability of being random associations and therefore 
this is the largest sample of statistically secure Spitzer counterparts 
to submm galaxies. The remaining two source have multiple pos- 
sible IRAC identifications, which would require somewhat deeper 
mid-IR or radio data to disentangle. 



The median redshift for our secure counterparts with opti- 
cal spectroscopic or photometric redshifts (21/35 sources) is 2.0, 
while the median redshift for all counterparts, including those with 
only IR photometric redshifts (33/35 sources), is 2.2. With conser- 
vative assumptions we place a limit of < 14 per cent on the fraction 
of 850 fim selected galaxies lying at z > 4. This is a lower frac- 
tion than is suggested for galaxies selected at ~ 1 mm (Eales et 
al. 2003). 

Having a large identified submm sample has allowed us to in- 
vestigate properties of the bright (Ssso > 2 mjy) submm popula- 
tion. In most cases, the submm flux is clearly dominated by a single 
24 fim source (20-700 /ily) with a wide range of 24 fim to 850 
flux ratios as a function of redshift, suggesting the presence spectral 
features and extinction in the mid-IR. 

We find significant scatter in the SEDs of submm galaxies, 
due to both a range of luminosities and SED shapes (tempera- 
tures). A composite rest-frame SED shows that the submm sources 
have SEDs which peak at longer wavelengths than those of local 
ULIRGs. Using a basic greybody model, 850 fim selected galax- 
ies appear to be cooler (~ 30 K) than local ULIRGs of the same 
luminosity. This may be telling us that the emission coming from 
high redshift submm galaxies is more extended than that in local 
ULIRGs, although it could also be due to selection effects both lo- 
cally and at high redshift. 

We also investigate the IR luminosities of the submm sources 
in GOODS-N, as estimated through their mid-IR, submm and radio 
flux densities. When all three measures are combined we find val- 
ues for LIR which are consistent with the values derived from the 
mid-IR and radio data separately. Using templates based on local 
ULIRGs with a fixed luminosity-temperature relation, the submm 
flux will typically overestimate the luminosity. However, when the 
temperature and luminosity are allowed to vary, 850 ^m, 24 
and 1.4 GHz flux density together give a good fit to the luminosity. 
The median IR luminosity (using radio, submm and infrared flux) 
for our sample is L (8 — 1000 fim) = 6.0 x 10 12 Lq, which cor- 
responds to a SFR of 1100 MQyr -1 . This is consistent with the 
values derived in Chapman et al. (2005) for 73 spectroscopically 
identified radio-detected submm galaxies using only the radio and 
submm flux densities. 

We compare the infrared properties of these submm galaxies 
to a large sample of BzK galaxies in GOODS-N (Daddi et al. 2005). 
The median 24 fim to 850 fim flux density ratio of submm galaxies 
is a factor of ~ 4 lower than that of their high redshift neighbours, 
the BzK galaxies, indicating that these galaxy populations are char- 
acterized by different SED shapes. We may understand this differ- 
ence in terms of the evolutionary scenario advocated by Vega et 
al. (2005), in which the submm galaxies, with cooler temperatures, 
may be an earlier phase in the star formation than the BzK galaxies. 
But a simple comparison of surface densities and luminosities also 
shows that the submm galaxies are more extreme star-formers than 
the average BzK galaxy. 

The IR luminosity (derived from all three wavelength flux 
densities) correlates with both the submm flux density and the red- 
shift. While selection effects will affect this to some level, it is hard 
to understand this correlation without invoking strong evolution of 
ULIRGs with redshift. 

At shorter wavelengths, the mid-IR SED shows evidence for 
the presence of the stellar 1.6 fim 'bump'. The Spitzer photometry 
alone can be used to constrain the redshift with reasonable accu- 
racy, a(Az/(l + z)) = 0.07. Roughly 2/3 of the secure submm 
counterparts satisfy the criteria Ss.$ > 6 and Ss.s/Ss.e > 1.3, 
which has a low probability of occurring at random within the 
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8arcsec search radius in GOODS-N. Hence a simple cut on this 
IRAC colour-magnitude plane is an efficient way of finding coun- 
terparts to submm galaxies. 

Several sources in our sample show a power-law shape at mid- 
IR wavelengths, which may indicate the presence of a powerful 
AGN. However, when combined with deep X-ray imaging, we find 
that only 1/21 (5 per cent) of secure sources in our sample are 
likely to be dominated by an AGN in the mid-IR emission. For the 
bulk of the sources the IR luminosity appears to be dominated by 
star formation rather than nuclear activity. 

Spitzer has opened up a new wavelength range through which 
we can study submm galaxies. Not only has it allowed us to identify 
many more submm counterparts, but we can begin to explore some 
of their physical properties. In the future, deep IRS spectroscopy 
will provide a more detailed look at the power behind the intense 
star formation in submm galaxies and disentangle many of the is- 
sues discussed in this paper. 
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APPENDIX A: NOTES ON INDIVIDUAL SOURCES 

Here we summarize the identifications for the full sample of 40 
submm sources in GOODS-North, which are found at > 3.5<r in the 
SCUBA super-map. In Section|2| we discussed the flux-deboosting 
of our sample. The following Paper III sources are not included in 
the analysis of this paper due to the fact that they are very noisy 
and have a non-negligible probability (> 5 per cent) of having es- 
sentially zero flux density: GN27, GN29, GN33, GN36 and GN38. 
However, for completeness (and since it is likely that many of them 
are real sources), we also discuss the possible counterparts for these 
5 less secure sources here. Our procedure for counterpart identifi- 
cation is outlined in Section|4| 

Using the Spitzer data and new radio reduction, we are now 
able to find unique counterparts for 33/35 submm sources, 21 of 
these are secure counterparts (P < 0.05) and 12 are tentative coun- 
terparts (0.05 < P <0.20). The remaining two sources, GN08 and 
GN35, have multiple IRAC counterparts within the search region 
which are equally likely and therefore we are unable to assign a 
unique identification. All of the submm sources with previous VLA 
1.4 GHz detections (Richards 2000) were confirmed in the new re- 
duction (Morrison et al.,in preparation) and, in addition, there are 
four sources from the Paper III radio-detected list which now have 
a second radio source within the search radius, with three of these 
being close enough to be regarded as interacting systems. Of the 
12 optical counterparts suggested in Paper III for radio-undetected 
sources using the optical data, 7 are confirmed as being correct 
using the Spitzer data and the new radio images. These new data 
favour different counterparts for four of these sources, although the 
Paper III counterparts are still feasible in some cases and may make 
a sub-dominant contribution to the submm emission in others. The 
last source has several possible counterparts, including the Paper 
III counterpart. Six of these previously radio-undetected submm 
sources are now detected above 3cr in the 1.4 GHz radio image. 

Al Secure counterparts 

For each of the 21 sources described in this section, the probability 
that the counterpart is associated with the submm source at random, 
based on having a > 3<r radio detection, MIPS 24 /im detection, 
and/or red IRAC colours is < 5 per cent. We now discuss each 
of these sources in turn and refer to Fig. lAll which shows postage 
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stamp images for each source at optical, IRAC, MIPS and radio 
wavelengths. 

GN03: There is no MIPS/IRAC detection of the Paper III iden- 
tification. However, there is a relatively bright MIPS source, with 
a 4.3a radio detection, within the search radius. There is an IRAC 
detection of this source, although there is no ACS detection. 

GN04: There are two radio sources within the search ra- 
dius, separated by 2.5 arcsec, with one of them twice as bright as 
the other. Both have an IRAC and a hard X-ray counterpart, but 
the fainter radio source lacks an optical counterpart. There is one 
24 fim source located between the two, but it is difficult to separate 
this into two sources, given the 24 fim resolution. The original VLA 
source has a spectroscopic redshift of 2.6 (Chapman et al. 2005). At 
this redshift, the separation between the two radio sources is only 
~ 20 kpc. The IRAC flux densities of the two sources are compara- 
ble. While the ACS image does not show any interaction between 
them, the morphology of the brighter radio source is disturbed and 
has multiple components (Paper III). We conclude that the submm 
emission is coming from two interacting galaxies at the same red- 
shift. This submm source also has a close companion (25 arcsec 
away), GN04.2 discussed in Section A3. However, the counterparts 
of each seem to imply that they are at different redshifts and not part 
of the same system. 

GN05: The only MIPS source in the search radius is coin- 
cident with the Paper III identification and has a new 3.9a radio 
detection. 

GN06: There is one likely counterpart which is detected in 
radio, IRAC, MIPS, and ACS, making this an unambiguous identi- 
fication. 

GN07: There are two radio sources within the search radius, 
separated by 2.5 arcsec. One of these radio sources is three times 
brighter than the other. Both have an ACS, IRAC and MIPS coun- 
terpart. Interestingly the fainter radio source is brighter in the op- 
tical and mid-IR. The IRAC fluxes for the two sources are almost 
identical at all four wavelengths, and the shape of the IR SED is 
consistent with both sources being at the same redshift. Chapman 
et al. (2004a) report a spectroscopic redshift of 1.99 for the brigher 
optical, fainter radio source. At this redshift 2.5 arcseconds corre- 
sponds to about 20 kpc, so it is likely that this is an interacting or 
multi-component system. 

GN10: The Paper III identification has an IRAC detection, 
however there is a new 4.5ct radio detection 2.5 arcseconds from 
this system. The new radio source has a faint IRAC detection but 
no optical counterpart. There is very faint MIPS emission which 
extends between these two sources. The Paper III counterpart is at 
z — 1.344 and therefore it should be fairly bright at 24 /im if it 
is a submm source. Since it is not, we conclude that the new radio 
source is the most likely submm counterpart, but note that the Paper 
III counterpart and this new radio source could be associated. The 
optical counterpart from Paper III is classified as being asymmetric 
which is an indication of a merging, or disturbed, system. 

GNU: In Papers II and III this submm source was identified 
with a WSRT radio source (Garrett et al. 2000). However, the po- 
sitional uncertaintly of the WSRT source is ~ 3 arcseconds. There 
is no MIPS/IRAC detection of the WSRT source, but on the other 
hand there is a new radio source with a MIPS/IRAC counterpart 
almost 5 arcseconds from the WSRT position which is undetected 
in the optical. By fitting the 850/24 fim flux to a suite of templates, 
we find that neither of the two other MIPS sources in the search 
radius are likely to be the source of the submm emission, and be- 
cause the radio source has a low probability of being there at ran- 
dom, we assign the new radio source as the counterpart. This is the 



only radio-detected Paper III counterpart which has changed. This 
submm detection is also in the Chapman et al. (2005) catalogue, 
although the suggested identification there is different, and there is 
no detection of their counterpart in the new radio image. At a red- 
shift of 0.5, the source in the Chapman et al. (2005) catalogue is 
not likely to correspond to the submm source, given the lack of a 
bright VLA radio detection and low 24 fim flux. 

GN12: There is one likely counterpart, which is detected in 
radio, IRAC, MIPS and ACS, making this an unambiguous iden- 
tification. The MIPS image in this region is crowded and blend- 
ing is an issue, therefore there is a large systematic uncertainty in 
the 24 fim flux. We found that the addition of the IRAC photom- 
etry improved the photometric redshift estimate significantly, and 
therefore we use the new redshift in this paper (see Table A2). 
In the ongoing campaign to obtain more redshifts in GOODS-N, 
we have targeted several of the submm counterparts with Keck. 
However, due to poor weather and the faintness of the sources in 
the optical, many of the spectra have come up blank. The LRIS 
slit for GN12 happened to pick up another 24 fim galaxy ~ 8 
arcsec away from the SCUBA position. We detected strong Ly- 
man alpha emission at z = 2.006 from this serendipitous source 
(12 h 36 m 46.73 s , +62° 14'46.3"). While this is near the peak of the 
redshift range for submm counterparts, it is elimimated as a pos- 
sible counterpart because: 1) it is too far from the SCUBA po- 
sition; and, 2) it is bright in the ultraviolet, which is not a char- 
acteristic of submm counterparts. In fact, the submm observations 
of ultraviolet bright galaxies such as LBGs have shown that while 
both populations are star forming galaxies at similar redshifts, the 
bright submm and LBG populations overlap very little (Chapman 
et al. 2000; Webb et al. 2003). However, there is evidence for strong 
clustering between bright submm and LBG galaxies which might 
indicate that they lie at the same redshifts (Webb et al. 2003; Borys 
et al. 2004a). Therefore, while this new z — 2.006 24 fim galaxy 
is not likely to be producing the submm flux, it may be part of the 
same high redshift system. 

GN13: The counterpart to this source is detected in radio, 
IRAC, MIPS, and ACS. While there are two additional MIPS 
sources within the search radius capable of producing the 850 fim 
emission, neither of them have a radio detection, which might be 
an indicator of the correct identification. This submm counterpart 
has the lowest redshift in our sample, with z spoc = 0.475. The 
counterpart to this submm source is a supernova host (Chary et 
al. 2005), which is the highest IR luminosity SN host in GOODS- 
N at z < 0.5, and therefore it is not surprising that it is the only one 
we detect with SCUBA. At this low 850 fim flux level, it is reason- 
able for the submm emission to be associated with the supernova 
host galaxy (see Farrah et al. 2004). 

GN14: This source is also known as HDF850.1, and is dis- 
cussed in detail in several papers (see Dunlop et al. 2004, and ref- 
erences therein). There is a low redshift (z spcc = 0.300), radio 
detected galaxy 6 arcseconds to the south-west of the submm posi- 
tion with a P value of 0.04. While this would qualify as a secure 
counterpert under our criteria, there is additional evidence which 
suggests that the counterpart for GN14 is much more complicated. 
Due to a detection with IRAM and MERLIN, the counterpart to 
GN14 is thought to lie behind a foreground elliptical galaxy very 
near the submm position. In the new reduction of the VLA data, we 
find a 2.7ct peak at the position of the Dunlop et al. (2004) counter- 
part. The position of the elliptical and the submm counterpart are 
less than an arcsecond apart. There is an IRAC and MIPS detec- 
tion of this system, but given the resolution of the Spitzer data, it is 
difficult to determine the contributing flux from the submm coun- 
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terpart and elliptical separately. There is also an additional com- 
plication from the foreground elliptical galaxy possibly lensing the 
submm system. Nevertheless, the probability of randomly finding 
a red IRAC source which is also detected at 24 /im this close to 
the submm position is < 0.05. The faintness of the MIPS flux is 
consistent with the photometric redshift for the submm source be- 
ing ~4 (Dunlop et al. 2004). We list this source in Table lATI and 
Table lA2l and consider it a secure identification. However, since we 
can only put upper limits on the IRAC and MIPS fluxes, we leave 
it out of the figures and analysis in this paper. Given that the Pois- 
son statistics favor a different counterpart, GN14 is an example of 
a case where the strict statistical process fails. 

GN16: There is one likely counterpart which is detected in 
radio, IRAC, MIPS and ACS, making this an unambiguous identi- 
fication. 

GN17: There are two radio sources within the search radius, 
separated by 5.5 arcsec. One of them is three times brighter then 
the other. Both sources have an IRAC/MIPS counterpart, although 
the fluxes of the IRAC and MIPS counterparts are different by fac- 
tors of 5 and 10, respectively. The brighter radio source is much 
fainter optically and has a photometric redshift of 1.72 (Paper III) 
which would indicate a separation of 45 kpc if both sources are at 
the same redshift. This is perhaps a little too large to be considered 
one system and, in fact, the fainter radio source has a photometric 
redshift of 1.2. By fitting the MIPS flux to a range of templates we 
find that the brighter radio source is more likely to be the dominant 
submm emitter, and therefore assign this as the counterpart. Cohen 
et al. (2000) list a spectroscopic redshift for this galaxy of 0.884, 
but given its optical and IR colours, such a redshift seems very un- 
likely. This counterpart satisfies the 'BzK' criteria for classification 
as a massive galaxy at z > 1.4 (see Daddi et al. 2005 for more 
details), and therefore we use our photometric redshift rather than 
the apparently discrepant spectroscopic redshift. 

GN18: There is only one MIPS source in the search radius and 
it has an IRAC counterpart but no optical detection. There is also a 
new 3.7a radio detection of this source. 

GN19: There are two equally bright radio sources within the 
search radius, separated by 3.0 arcsec. Both have an IRAC/MIPS 
counterpart, but the brighter MIPS source lacks an associated op- 
tical source. Nevertheless, both of these objects are confirmed to 
have a spectroscopic redshift of 2.484 (Chapman et al. 2005). At 
this redshift, the separation is 24 kpc and therefore we conclude 
that the submm emission is coming from two interacting galaxies 
at the same redshift. 

GN20: The counterpart to this submm source is now unam- 
biguously known, due to detections with both IRAM and the SMA 
(Pope et al. in preparation; Iono et al. 2006). The new radio reduc- 
tion has also revealed a source at 1.4 GHz which looks extended. 
This counterpart is detected in MIPS, IRAC and ACS. The source 
is a _B-dropout, which implies that it is at high redshift, and we are 
currently obtaining optical spectroscopy to confirm this. 

GN20.2: The counterpart to this source is detected in radio, 
IRAC, MIPS and ACS. It is the companion to the bright submm 
source GN20 (18 arcsec away) and is also a B-dropout galaxy. 
Note that GN20.2 was first published in Chapman et al. (2001) 
when they targeted optically faint radio sources. However, so far 
it has eluded a spectroscopic identification (S. Chapman, private 
communication). 

GN22: There is one likely counterpart, which is detected in 
radio, IRAC, MIPS and ACS, making this an unambiguous identi- 
fication. The other MIPS source within the search radius which is 
coincident with an edge-on galaxy (see Fig. IA1I must contribute 



less than half of the submm emission since it is at lower redshift 
(z spcc = 0.6393). 

GN25: There is one likely counterpart, which is detected in 
radio, IRAC, MIPS and ACS, making this an unambiguous identi- 
fication. Note that the radio emission is extended (see Fig. lAlt . 

GN26: There is a bright radio, MIPS, IRAC and optical coun- 
terpart with a spectroscopic redshift of 1.219. Given the high flux 
density at radio and mid-IR wavelengths, this is considered a se- 
cure counterpart even though it has a relatively large offset from 
the submm position. 

GN30: This source was identified with a faint 8.5 GHz radio 
source in Paper II, although there was no detection at 1 .4 GHz in the 
Richards (2000) catalogue. In Paper III we found an optical coun- 
terpart at the 8.5 GHz position. There is an IRAC/MIPS detection 
of this galaxy, and in addition there is a 3a peak in the new 1.4 GHz 
radio image. This is a very blue, high surface brightness 'chain 
galaxy', more similar to typical LBGs than most submm counter- 
parts. However, given the low 850 /im flux, it is perfectly reason- 
able that this is the counterpart. There is another MIPS source in 
the search radius, although it is at lower redshift (z spac — 0.557). 
Given the radio and mid-IR detection, it qualifies as a secure coun- 
terpart. 

GN31: There is no MIPS or IRAC detection of the Paper III 
counterpart, but there is one relatively bright MIPS source in the 
search radius, which has a 3<r peak in the new radio image. This 
MIPS source is detected with both IRAC and ACS, and it is one of 
the lowest redshift counterparts, z spcc = 0.935. 

A2 Tentative counterparts 

The following 12 sources have likely radio, 24 fj,m, and/or IRAC 
counterparts but, because the probability of random association is 
above 0.05, the identification is listed as tentative. 

GN01: There is only one MIPS source in the search radius. It 
is coincident with the Paper III optical identification and has a 3.5a 
detection in the new radio map. 

GN02: There is a new 5.6a radio source just beyond the search 
radius. The new radio source has a bright MIPS, IRAC and ACS 
counterpart. There are no other MIPS sources within the search 
radius, and we therefore assign this radio source as the counterpart. 
This is the only proposed counterpart which is beyond the search 
radius, but since the submm noise is high, it is not unreasonable for 
it to have a greater positional offset. 

GN04.2: There is only one MIPS source in the search radius 
and it has an IRAC and ACS detection and a new 3a radio detec- 
tion. This submm source is the companion to GN04. However, the 
redshift estimates of the counterparts seem to imply that they are 
independent systems. 

GN09: There are two MIPS/IRAC sources within the search 
radius. By fitting the 850/24 flux to a suite of templates, we 
find that the second source is fainter at 24 /im and is not capable of 
contributing more than 10 per cent of the submm flux. In addition 
the brighter 24 /im source has red IRAC colours, which means that 
it has a low probability of being within the search radius at random. 
We therefore assign the brighter MIPS source as the counterpart. 
This counterpart is not detected in the ACS or radio images, and is 
therefore likely to be at high redshift. 

GN15: There is a MIPS/IRAC detection of the Paper III coun- 
terpart. This source is bright at 24 /im and has red IRAC colours. 
The VLA source to the North, discussed in Paper III, is detected 
in the new radio image and has a P value of 0.05. Although there 
appears to be extended MIPS emission close to the radio source, 
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Figure Al. 'Postage stamp' images of submm counterparts, from left to right: ACS zgso; IRAC 3.6 /im; MIPS 24 fim; and VLA 1.4 GHz. All images are 
20 X 20 arcsec and the black circle is centred on the SCUBA position with a 8 arcsec radius. The smaller square, or circle, indicates the secure, or tentative, 
counterpart, respectively. 
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it is very faint, and moreover, there is no IRAC or ACS galaxy as- 
sociated with it, although there are several nearby. All of the other 
submm galaxies which are radio detected are bright at 24 /im and 
therefore the mid-IR faint, radio source is not typical of submm 
galaxy counterparts. Since the 24 /im flux of the Paper III counter- 
part is so high there is a lower probability that it is a random asso- 
ciation. This is a rare case where a radio detection is not considered 
the most likely counterpart. 

GN21: There is only one MIPS source in this region and it has 
a 4.7a radio detection and a red IRAC counterpart. This source is 
undetected in the optical and the IR photometric redshift puts it at 
3.3. There are no other likely counterparts in the search radius. 

GN23: This source has a radio, IRAC, and MIPS detection. 
However, it is only a faint smudge in the ACS image, and therefore 
we are not able to obtain an accurate optical photometric redshift, 
but we do have an estimate of the redshift from the Spitzer photom- 
etry (see Section lo31 and Table A2). 

GN24: There is only one MIPS source in the search radius 




continued... 



and it is coincident with the Paper III counterpart. There is no radio 
detection but it is coincident with a red IRAC source. 

GN28: This source is unique in that it is within 8 arcseconds 
of a massive radio jet (Borys et al. 2002). Within the search ra- 
dius, there is faint MIPS emission, elongated in the NE-SW di- 
rection, which connects several IRAC sources. One of these IRAC 
sources is 0.6 arcseconds from the end of a radio jet and so the 
usual assumption that radio flux is indicative of the submm coun- 
terpart is more difficult here. The ~ 15 arcsec long radio feature 
is the most prominent in the entire GOODS -N field (Richards et 
al. 1998). There is a disturbed optical detection coincident with 
the end of the radio jet, which is confirmed to be at a similar red- 
shift (2 spec = 1.020, Phillips et al. 1997) to that of the central 
AGN (z spcc = 1.013, Richards et al. 1998). The fact that the 
submm emission seems coincident with the end of the jet is sug- 
gestive of jet-induced star formation. This is not the only 850 /im 
selected source which lies near a radio jet. LE850.07 (also known 
as LH1200.096, Ivison et al. 2002; 2005) is a submm source as- 
sociated with a radio lobe, however spectroscopy of the lobe and 
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the core fail to put them at the same redshift and there are other 
candidate counterparts in the field. In both of these cases, it will 
be necessary to better localize the submm emission before draw- 
ing firm conclusions. There is also another IRAC source with faint 
MIPS emission which is an optically-invisible X-ray source in the 
region. While there are several possible counterparts, the fact that 
other submm sources have been seen close to radio jets, which has 
a low probability of happening at random, convinces us to assign 
the MIPS/IRAC source at the end of the radio jet as a tentative 
counterpart. 

GN32: The Paper III optical counterpart has a MIPS and IRAC 
detection, but it is not detected in new radio image. The 24 ^im flux 
is fairly high and it has red IRAC colours. The only other MIPS 
source in the search radius is radio-detected and known to be at 
lower redshift (z sp cc = 0.761). This radio/MIPS source has a P 
value of 0.06, therefore appears statistically likely to be the submm 
counterpart. However, based on reasonable SEDs which are consis- 
tent with the rest of our sample, this source is not likely to dominate 
the submm flux, although of course it may contribute. Because of 
this we maintain that the Paper III source is the major submm emit- 
ter. 

GN34: The Paper III counterpart has a MIPS and IRAC de- 
tection, however there are also two other MIPS sources within 
the search radius, both of which are capable of contributing to 
the submm flux. Based on Poisson statistics, the red IRAC/MIPS 
source nearest the submm centre has, by far, the lowest probability 
of random association and therefore we assign it as the tentative 
counterpart. 

GN37: The Paper III counterpart has a MIPS and IRAC detec- 
tion, but it not detected in the new radio image. The MIPS flux is 
fairly low but it has red IRAC colours. 

A3 Multiple counterparts 

The following two submm sources have multiple possible counter- 
part within the search radius which are equally likely. Since there is 
not one single counterpart which sticks out as having a significantly 
low probability of random association, we are unable to assign a 
unique counterpart for these sources. The submm emission from 
these sources could be coming from one of the multiple possible 
counterparts or from several of them. 

GN08: There is an IRAC detection of the Paper III counter- 
part, but no MIPS or radio detection. However, the lack of 24 ^im 
detection may be because it is in the wings of much brighter MIPS 
sources, which are outside the search radius and therefore we can 
only obtain an upper limit on the 24 pirn flux density. There are two 
MIPS detected sources within the search radius but neither of these 
can contribute more than 10 percent of the submm flux, since they 
are at low redshift and faint at 24 fim. Since the probability of ran- 
dom association is not significantly low for either of these sources, 
we are unable to assign a unique identification. It will be possible 
to make further progress on this by localizing the submm emission 
using IRAM or SMA interferometry. 

GN35: There are no MIPS sources in the search radius, al- 
though there are three faint IRAC/ACS sources. Two of these can 
be eliminated as being likely counterparts, since their optical pho- 
tometric redshifts put them at z < 0.1. The third IRAC source is 
detected in the ACS image and has a photometric redshift of 2.20. 
Given that all other submm sources in GOODS-N are detected in 
IRAC, there is no reason to believe that this one is not. However, 
since the probility of random association is not significantly low, 
we are unable to assign a unique identification. 



A4 Possible spurious sources 

The following sources are no longer in our sample, due to the 
fact that each of them individually has a non-negligible probabil- 
ity (> 5 per cent) of being at essentially zero flux density once we 
apply flux-deboosting (see Coppin et al. 2005). However, we dis- 
cuss their possible counterparts here for completeness and because 
the chances are high of several of them being real. 

GN27: There is only one MIPS source in the search radius, 
and it has an IRAC and ACS counterpart. This is the most likely 
identification. 

GN29: There is only one MIPS source in the search radius 
and it has an IRAC and ACS counterpart. This is the most likely 
identification. 

GN33: There are three possible IRAC sources and no MIPS 
sources in the search radius, therefore we cannot assign a unique 
identification. 

GN36: There are no MIPS sources in the search radius. There 
are three IRAC/ACS sources, but nothing that stands out and there- 
fore we cannot assign a unique counterpart. 

GN38: This source is near the edge of the original GOODS- 
N region, therefore there is no deep ACS coverage here, although 
there are radio and Spitzer data. There are no Richards (2000) VLA 
sources within the search radius, but there are three sources in the 
new radio reduction. The furthest radio source is coincident with a 
Zspcc = 0.1136 galaxy and therefore is not likely to be the submm 
counterpart. The two closer radio sources have similar radio flux, 
are separated by 4.5 arcseconds, and each of them is associated with 
an IRAC and MIPS source of comparable flux. Using just the in- 
frared, submm and radio information, these two sources have very 
similar SEDs, which may indicate that they are at the same red- 
shift. If this is true then these sources are both contributing to the 
submm flux (approximately equally in this case), and are possibly 
an interacting system at high redshift. Given the three other submm 
sources in our sample with two radio counterparts, we consider this 
to be the most likely scenario and consider both radio sources as 
the correct identification. Since they are radio-bright and close to 
the submm position, this system would classify as a secure coun- 
terpart. 
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Table Al. Submm sources in GOODS-North: submm, and counterpart, positions, and submm flux densities. The second column lists the full submm name 
which indicates the SCUBA position. The third and forth columns list the position of the IRAC counterpart and have been adjusted for the known offset to 
the VLA image (Dec(VLA) - Dec(GOODS-N) = - 0.38 arcsec). The last two columns give the raw and corrected 850 ^m flux density after applying the 
flux deboosting algorithm of Coppin et al. (2005). We have not applied a correction to sources with high signal-to-noise ratio (> 6) or a very low noise level 
(< 0.7mJy), since the correction to these sources is negligible compared to our calibration accuracy. After rejecting the five sources with very high noise 
levels, all submm sources have a likely unique counterpart, except for GN08 and GN35, where there are multiple possible counterparts. 



SMM ID 


.S IV1 IV1 Name 


IRAC position 


Raw Ssso/im 


.... . . . , 1 o 

Corrected b &50ilrn 






RA 


Dec. 


(mjy) 


/.„ T,,\ 

(mJy) 


GN01 


SMMJ123606.7+621556 


12-36-06 70 


62-15-50 43 


7.3 ±1.5 


6.2 ±1.6 


GN02 


SMMJ123607.7+621147 


12-36-08 81 


62-11-43 57 


16.2 ±4.1 


12.1 ±4.3 


GN03 


SMMJ123608.9+621253 


12:36:08.65 


62:12:50.81 


16.8 ±4.0 


12.8 ±4.2 


GN04 a 


SMMJ123616.6+621520 


12:36:16.11 


62:15:13.53 


5.1 ± 1.0 


4.9 ±0.7 






12:36:15.82 


62:15:15.34 






GN05 


SMMJ123618. 8+621008 


12:36:19.13 


62:10:04.32 


6.7 ± 1.6 


5.2 ±1.8 


GN06 


SMMJ123618.7+621553 


12:36:18.33 


62:15:50.40 


7.5 ±0.9 


7.5±0.9 b 


GN07 a 


SMMJ123621.3+621711 


12:36:21.27 


62:17:08.16 


8.9 ±1.5 


8.9±1.5 b 






12:36:20.98 


62:17:09.55 






GN08 


SMMJ123622.2+621256 






12.5 ±2.7 


10.0 ±3.0 


GN09 


SMMJ123622.6+621617 


12:36:22.07 


62:16:15.81 


8.9 ± 1.0 


8.9 ± 1.0 b 


GN10 


SMMJ123633.8+621408 


12:36:33.40 


62:14:08.72 


11.3 ± 1.6 


11.3 ± 1.6 b 


GNU 


SMMJ123637.2+621156 


12-36-37 51 


62-1 1-56 52 


7.0 ±0.9 


7.0 ± 0.9 b 


GN12 


SMMJ123645.8+621450 


1 l-'Kf.-Afx 07 

1Z,.jU.HU.U / 


(S?- 14-48 7(S 


8.6 ± 1.4 


8.6±1.4 b 


GN13 


SMMJ123650.5+621317 


1 7-^fv4Q 77 


AO-i Q. 1 7 on 


1.9 ±0.4 


1.9±0.4 b 


GN14 


SMM J 123652.2+62 1226 


lZ,.J\J.JZ,.\JO 


fS7- 1 7-76 71 


5.9 ±0.3 


5.9±0.3 b 


GN15 


SMM J 123656.5+62 1202 


IZ.JU.JJ.Oi 




3.7 ±0.4 


3.7±0.4 b 


GN16 


SMM J 123700.4+6209 11 


i 76 

1Z..J/ .\J\J.Z.\J 


(S7-OQ-OQ 77 


9.0 ±2.1 


6.9 ±2.5 


GN17 


SMMJ123701.2+621147 


17-^7-01 SO 

1Z.J/.U1.J7 


fV7- 1 1 -Af\ 7S 
OZ. 1 1 .H-O.Z J 


3.9 ±0.7 


3.9±0.7 b 


GN18 


SMMJ123703.0+621302 


1Z.J /.UZ.JJ 


A7- 1 ^-(17 77 
DZ. 1 J.UZ.ZZ 


3.2 ±0.6 


3.2±0.6 b 


GN19 a 


SMMJ123707.7+621411 


iz.j/.u/.iy 


A7- 1 A-(\l 07 
OZ. l^f.U / .7 1 


10.7 ±2.7 


8.0 ±3.1 






1 7 -^7 -07 


A7- 1 A-no. 
oz. I'+.uy. ju 






GN20 


SMMJ12371 1.7+622212 


1 7* IT- 1 1 SC 
1Z. J /.I 1.85 


^0-09- 10 11 
OZ.ZZ. 1Z. 1 1 


20.3 ±2.1 


20.3 ± 2.1 b 


GN21 


SMMJ123713.3+621202 


12:37:14.06 


62:11:56.75 


5.7 ± 1.2 


4.9 ± 1.2 


GN22 


SMMJ123607.3+621020 


12:36:06.84 


62:10:21.36 


14.4 ±3.9 


10.5 ±4.3 


GN23 


SMMJ123608.4+621429 


12:36:08.60 


62:14:35.30 


7.0 ± 1.9 


4.9 ±2.3 


GN24 


SMM J 1 236 1 2.4+62 1217 


12:36:12.00 


62:12:21.99 


13.7 ± 3.6 


10.0 ± 4.1 


GN25 


SMMJ123628. 7+621047 


12:36:29.12 


62:10:45.92 


4.6 ± 1.3 


3.2 ± 1.4 


GN26 


SMMJ123635.5+621238 


12:36:34.51 


62:12:40.93 


3.0 ± 0.8 


2.2 ± 0.8 


GN28 


SMMJ123645. 0+621 147 


12:36:44.51 


62:11:41.92 


1.7 ± 0.4 


1.7 ± 0.4 b 


GN30 


SMMJ123652.7+621353 


12:36:52.75 


62:13:54.59 


1.8 ± 0.5 


1.8 ± 0.5 b 


GN31 


SMMJ123653. 1+621 120 


12:36:53.22 


62:11:16.69 


2.8 ± 0.8 


2.1 ± 0.6 


GN32 


SMMJ123659. 1+621453 


12:36:58.74 


62:14:58.92 


5.3 ± 1.4 


3.8 ±1.6 


GN34 


SMMJ123706.5+622112 


12:37:06.22 


62:21:11.57 


5.6 ± 1.6 


3.8 ±1.9 


GN35 


SMMJ123730.8+621056 






14.3 ±3.9 


10.4 ±4.3 


GN37 


SMMJ123739. 1+621736 


12:37:38.26 


62:17:36.38 


6.8 ± 1.9 


4.7 ±2.3 


GN04.2 


SMMJ123619.2+621459 


12:36:18.67 


62:15:03.09 


3.6 ± 1.0 


2.7 ± 1.0 


GN20.2 


SMMJ123709.5+622206 


12:37:08.77 


62:22:01.78 


11.7 ± 2.2 


9.9 ±2.3 


GN27 


SMMJ123636.9+620659 


12:36:36.00 


62:07:00.51 


24.0 ±6.1 


17.6 ± 5.8 C 


GN29 


SMMJ123648.3+621841 


12:36:48.17 


62:18:42.93 


20.4 ±5.7 


14.7 ± 5. 6 C 


GN33 


SMMJ 123706.9+62 1850 






21.8 ±5.8 


15.9 ± 5.6 C 


GN36 


SMMJ12373 1.0+621856 






24.8 ±7.0 


17.7±6.5 C 


GN38 a 


SMMJ123741.6+621226 


12:37:41.64 


62:12:23.61 


24.9 ± 6.5 


18.1 ±6.1 c 






12:37:41.16 


62:12:20.89 







a These submm sources appear to have two radio/IRAC counterparts contributing to the submm flux. 
b No correction is applied to these sources, since they are either high signal-to-noise or have very low noise levels. 
c These five sources are rejected from our sample based on the fact that they have a non-negligible probability of deboosting to zero flux 
density; however, we include them here since the probability is high that most of them are real. 
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Table A2. Multi-wavelength photometry of submm counterparts. There are 21 secure counterparts, and 12 tentative counterparts (see Section |3] for more 
details). We expect ~ 1 secure submm counterpart to be a random association from both the individual P values and Monte Carlo simulations. The tentative 
counterparts have a higher probability of being random associations. Within each of these two groups, sources are ordered according to their redshift, as 
measured from optical spectroscopy or photometry. Bold-faced redshifts are spectroscopic, otherwise they are photometric. Sources with no optical redshift 
information have tentative redshifts determined from IRAC colours listed in brackets (see Section l6"3l . Radio flux densities are listed if the detection is > 3<r, 
otherwise the 3<r upper limit from the local RMS is quoted. The second to last column lists the probability that the counterpart is a random association. The 
letter after the probability indicates which catalogue the random probabilities are based on: R (> 4cr radio), M (> 3cr MIPS 24 fim), R/M (> 3cr radio and 
> 3<t MIPS 24 fim), M/I (MIPS and red IRAC detection). See Section0foi more details on these probabilities. In the last column we list the IR luminosity 
determined from fitting the mid-IR, submm and radio flux densities to a suite of modified CE01 templates. LIR is defined as the integral under the best-fit 
rest-frame SED from 8-1000 fim. Note that we only list LIR if the source is detected at 1.4 GHz and 24 fim and has a reliable redshift estimate. 



SMM ID 


z 


'775 ma 8 


S 3.6 


S 4.5 


S5.8 


s 8.0 


S24 


Si. 4 


P 


LIR 






(AB) 




(MJy) 


(pJy) 




(M^y) 


(M^y) 




(10 12 L Q ) 


Secure submm counterparts 


GN13 C 


0.475 


21.6 


39.43 ± 1.24 


39.94 ± 0.99 


36.53 ± 1.49 


100.28 ± 1.74 


371.0 ± 10.4 


45.4 ± 5.4 


0.042 R/M 


0.13 


GN31 


0.935 


21.8 


60.64 ± 1.24 


44.83 ± 0.99 


39.91 ± 1.49 


34.64 ± 1.74 


367.0 ± 6.4 


16.3 ± 5.4 


0.034 M 


0.31 


GN25 C 


1.013 


22.8 


92.74 ± 1.24 


85.64 ± 0.99 


75.42 ± 1.49 


73.02 ± 1.74 


724.0 ± 12.0 


93.8 ± 12.9 


0.020 R 


1.1 


GN26 C 


1.219 


22.7 


63.11 ± 1.24 


72.06 ± 0.99 


59.55 ± 1.49 


75.32 ± 1.74 


446.0 ± 5.1 


194.3 ± 10.4 


0.027 R 


1.6 


GN30 


1.355 


22.7 


14.29 ± 0.89 


14.61 ± 0.79 


10.92 ± 1.28 


12.80 ± 1.48 


100.0 ± 5.9 


19.9 ± 5.4 


0.024 RIM 


0.36 


GN16 


1.68 


26.7 


17.04 ± 0.89 


25.40 ± 0.99 


32.64 ± 1.49 


26.25 ± 1.48 


267.0 ± 7.7 


361.9 ± 6.0 


0.00093 R 


5.2 


GN17 


1.72 


27.7 


53.87 ± 1.24 


70.13 ± 0.99 


71.54 ± 1.49 


57.95 ± 1.74 


710.0 ± 7.8 


110.0 ± 11.0 


0.0072 R 


5.5 


GN10 


(1.8) 


> 28 


1.21 ± 0.39 


1.96 ± 0.36 


2.72 ± 0.88 


5.11 ± 1.14 


30.7 ± 5.4 


33.1 ± 10.3 


0.027 R 




GN06 


1.865 


27.4 


14.54 ± 0.89 


19.41 ± 0.79 


26.94 ± 1.49 


20.07 ± 1.48 


330.0 ± 7.6 


178.9 ± 6.4 


0.0077 R 


7.2 


GN07 a 




27.8 


14.24 ± 0.89 


19.37 ± 0.79 


26.05 ± 1.49 


19.18 ± 1.48 


20.0 ± 6.2 


155.7 ± 6.3 


0.0055 R 


10.9 




1.998 


23.9 


14.82 ± 0.89 


18.67 ± 0.79 


25.10 ± 1.49 


18.58 ± 1.48 


346.9 ± 8.8 


55.10 ± 12.6 


0.016 R 




GN03 


(2.0) 


> 28 


5.29 ± 0.55 


8.27 ± 0.79 


12.69 ± 1.28 


14.56 ± 1.48 


172.0 ± 6.8 


29.1 ± 6.0 


0.022 RIM 




GN18 


(2.2) 


> 28 


2.53 ± 0.39 


4.12 ± 0.51 


6.66 ± 0.98 


7.63 ± 1.14 


79.9 ± 6.7 


21.2 ± 5.5 


0.050 RIM 




GNU 


(2.3) 


> 28 


9.71 ± 0.89 


13.89 ± 0.79 


19.56 ± 1.28 


19.39 ± 1.48 


145.0 ± 4.2 


25.6 ± 5.4 


0.016 R 




GN19 a 


2.484 


25.4 


17.97 ± 0.89 


23.30 ± 0.99 


31.76 ± 1.49 


28.95 ± 1.74 


34.7 ± 10.8 


39.0 ± 5.5 


0.037 R 


9.6 




2.484 


> 28 


9.83 ± 0.89 


14.51 ± 0.79 


22.23 ± 1.28 


22.51 ± 1.48 


242.7 ± 10.8 


38.0 ± 9.7 


0.0071 R 




GN22 


2.509 


24.6 


27.28 ± 1.24 


30.70 ± 0.99 


35.80 ± 1.49 


27.50 ± 1.74 


70.1 ± 5.2 


65.4 ± 11.8 


0.020 R 


3.6 


GN04 a 


2.578 


26.2 


12.32 ± 0.89 


18.07 ± 0.79 


29.49 ± 1.49 


43.36 ± 1.74 


302.8 ± 6.5 d 


58.0 ± 6.3 


0.049 R 


8.0 






> 28 


14.92 ± 0.89 


19.47 ± 0.79 


27.93 ± 1.49 


27.09 ± 1.74 




31.5 ± 6.3 


0.075 R 




GN05 


2.60 


24.9 


12.63 ± 0.89 


16.58 ± 0.79 


21.79 ± 1.28 


18.26 ± 1.48 


215.0 ± 6.0 


50.3 ± 15.4 


0.034 R/M 


6.5 


GN202 


2.83 


24.7 


3.87 ± 0.55 


3.89 ± 0.51 


6.06 ± 0.98 


9.36 ± 1.14 


30.2 ± 5.6 


180.7 ± 8.4 


0.028 R 


14.0 


GN20 


2.95 


24.4 


6.79 ± 0.55 


9.23 ± 0.79 


15.93 ± 1.28 


25.28 ± 1.48 


68.9 ± 4.8 


70.0 ± 16.3 


0.0092 R 


14.0 


GN12 C 


3.10 


26.2 


5.75 ± 0.55 


7.89 ± 0.79 


11.22 ± 1.28 


13.84 ± 1.48 


67.0 ± 9.0 


103.2 ± 5.5 


0.0056 R 


19.5 


GN14 b 


4.1 


> 24.6 


< 19.5 


< 17.6 


< 13.4 


< 16.2 


< 26.1 


< 15.9 


N/A 




Tentative submm counterparts 


GN04.2 


0.851 


22.6 


12.41 ± 0.89 


9.28 ± 0.79 


6.21 ± 0.98 


4.76 ± 1.14 


25.4 ± 5.2 


17.8 ± 5.9 


0.14 R/M 




GN2S 


1.02 


23.0 


9.24 ± 0.89 


0.54 ± 0.51 


5.68 ± 0.98 


4.35 ± 1.14 


20.9 ± 5.3 


90.3 ± 24.4 


0.080 R 


0.016 


GN02 C 


1.32 


23.8 


65.42 ± 1.24 


69.17 ± 0.99 


54.51 ± 1.49 


57.08 ± 1.74 


265.0 ± 6.2 


43.4 ± 6.1 


0.11 R 


1.4 


GN34 


1.36 


22.9 


18.4 ± 0.89 


18.7 ± 0.79 


14.0 ± 1.28 


18.3 ± 1.48 


81.7 ± 4.1 


< 23.6 


0.065 M 




GN32 


(2.1) 


27.8 


6.95 ± 0.55 


9.74 ± 0.79 


10.78 ± 1.28 


9.22 ± 1.14 


128.0 ± 6.2 


< 16.8 


0.15 M/I 




GN23 


(2.3) 


> 28.0 


6.45 ± 0.55 


9.51 ± 0.79 


13.44 ± 1.28 


18.31 ± 1.48 


55.9 ± 6.7 


34.9 ± 6.1 


0.10 R/M 




GN09 


(2.4) 


> 28 


3.97 ± 0.55 


6.60 ± 0.51 


10.56 ± 1.28 


14.18 ± 1.48 


49.8 ± 6.1 


< 19.5 


0.17 M/I 




GN01 


2.415 


23.3 


9.50 ± 0.89 


12.97 ± 0.79 


19.66 ± 1.28 


26.88 ± 1.74 


119.0 ± 6.0 


26.7 ± 12.1 


0.065 R/M 


4.9 


GN21 


(2.6) 


> 28 


3.95 ± 0.55 


5.64 ± 0.51 


9.29 ± 1.28 


10.05 ± 1.14 


46.3 ± 5.4 


33.5 ± 5.6 


0.11 R/M 




GN15 


2.743 


24.3 


15.68 ± 0.89 


18.74 ± 0.79 


22.04 ± 1.28 


18.95 ± 1.48 


200.0 ± 6.0 


< 16.2 


0.10 M/I 




GN24 


2.91 


24.7 


5.94 ± 0.55 


7.15 ± 0.51 


10.50 ± 1.28 


7.37 ± 1.14 


75.5 ± 6.2 


< 18.0 


0.12M/I 




GN37 


3.190 


23.1 


6.20 ± 0.55 


6.58 ± 0.51 


8.45 ± 0.98 


9.22 ± 1.14 


30.5 ± 4.6 


< 21.0 


0.12 M/I 





a These submm sources have two radio/IRAC counterparts very close together, both probably contributing to the submm flux, and therefore we list the multi-wavelength flux densities for both components. In all 

cases, both radio/IRAC sources are thought to be at the same redshift. 
^ D The Spitzer counterpart for this source (also known as HDF850. 1 ) is hidden behind a foreground elliptical galaxy, and therefore we quote the upper limits to the IRAC and MIPS flux densities. The 

photometric redshift of this source is from Dunlop et al. (2004). 
c These sources are also detected in 16 fim imaging of GOODS-N with the Spitzer IRS (Teplitz et al. 2005). 
d This MIPS source appears to be located between 2 IRAC sources. 
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